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ABSTRACT 


This thesis develops a decision aid to assist in assessing 
the cost effectiveness of upgrading a subsystem. The 
procedures developed in this thesis are to estimate the time 
of onset and the magnitude of the degradation of a subsystem 
and to estimate the best time to upgrade the subsystem. Two 
procedures are considered to estimate the time of onset of 
subsystem degradation and the magnitude of the degradation. 
One is maximum likelihood; the other is a Bayesian procedure. 
These estimates are then used in a cost model to estimate the 
cost of remaining with the current subsystem for the remaining 
planned lifetime of the system. A comparison of this cost with 
that of investing in the upgraded subsystem can be used to 
obtain a best time to invest in the upgraded subsystem. 
Procedures to assess the uncertainty of the cost advantage of 
upgrading the subsystem are also studied to give further 


information to the decision maker. 
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I. INTRODUCTION 


A. BACKGROUND 

A subsystem of an aircraft total system, for example an 
APS-80 antenna subsystem within the P-3B aircraft system, may 
tend to exhibit unfavorable reliability or maintenance cost 
properties, beginning at some random point in time. The 
evidence of such degradation of a subsystem suggests the 
possible economic and operational value of subsystem upgrade. 
The decision to upgrade a subsystem will, at least partially, 
be based on a comparison of the costs of remaining with the 
current subsystem over its remaining time horizon, and those 
of investing in the upgraded subsystem for use in the 
remaining time horizon. 

The cost of remaining with the current subsystem must be 
estimated using available data. It will depend on the estimate 
of the time of onset of subsystem degradation and the estimate 
of the magnitude and evolution of the degradation over time. 
The costs of investing in the upgraded subsystem are obtained 
from another source. A comparison of the cost (1) of remaining 
with the current subsystem with (2) investing in the upgraded 
subsystem can be used to obtain a "best" time to invest in the 


upgraded subsystem. This time can be beyond the system’s 


remaining planned life, in which case the upgrade is 


unadvisable. 


B. CURRENT METHOD AND MODEL 

The economic analysis program (ROI) in the Automated 
Management Indicator System (AMIS) is designed to compare an 
existing (current) subsystem with an improved (upgraded) 
subsystem. The program compares the projected cumulative cost 
of the existing subsystem to the projected cost of the 
improved subsystem. After computing these costs, the program 
computes the number of months until the cumulative cost of the 
existing subsystem is equal to the cumulative cost of the 
improved subsystem when the investment: cost and upgrade 
schedule are considered. If this break-even time exists, the 
time 1S printed on the screen. If this break-even point occurs 
before the end of the planned system life then the gains from 
upgrade are potentially available; if not then the upgrade is 
not likely to be worthwhile. 

Most of the input data for the economic analysis program 
comes from the NALDA data base. While this data base collects 
the measures of subsystem performance over a long period of 
time, the ROI program only uses an average of the data for the 
last 24 months to estimate the predicted future subsystem 
performance. It uses this average to compute the future cost 


of the current subsystem. 


Two features of the ROI model are considered in the 
development of a new model. One feature is that average values 
of measures of subsystem performance may be correct for local 
cost estimation, but may not represent the evolution of the 
subsystem degradation. The other feature is that the ROI 
program only computes the required time to cover the upgrade 
investment if the decision is made to upgrade the subsystem 
immediately. It may be better to wait to initiate subsystem 
upgrade in order to make more certain that an adverse trend 
exists. The model investigated in this thesis gives a best 


time to initiate subsystem upgrade based on estimated costs. 


C. NEW APPROACH 

The approach of this thesis is to consider the time series 
of a measure of subsystem performance and to estimate the time 
of onset of subsystem degradation and the magnitude of the 
degradation as time advances. These estimates are then used in 
a cost model to estimate the cost of remaining with the 
current subsystem for the remaining planned lifetime of the 
subsystem. We call this planned lifetime the time horizon. 

Two procedures are considered to estimate the time of 
onset of subsystem degradation and the magnitude of the 
degradation. Both are based on a simple change-point model 
that assumes that the degradation may begin at some time point 
and increase linearly thereafter. The change-point and the 


rate of degradation must be estimated from data; two 


procedures are used. One iS maximum likelihood; the other is 
a Bayesian procedure. Details of these procedures appear in 
Appendices B and C. A cost model is formulated to be 
compatible with that used in the AMIS ROI program. A 
description of the cost model appears in Appendix D. 
Procedures to assess the variability of the cost advantage of 
upgrading the subsystem are studied. The procedures used are 
described in the Appendix E. 

A decision aid using the methodologies is programmed in 
TURBO PASCAL. The resulting program is called UPGRADE.PAS. The 
program’s source code and user documentation are attached in 
Appendix H. The program can simulate data and provide 
graphical output. 

The methodologies are used to analyze simulated data and 
data from a radar transmitter on the F-14A. The data were 
supplied by C. Wrestler of NAVAIR (419) and come from the 
NALDA data base. These data appear in Appendix F. Results from 
the analysis appear in Appendix G. 

A brief discussion of the models and procedures appears in 
the next chapter. Information concerning the input data for 
this program appears in Chapter III. The results of analyses 
conducted using the PASCAL program will be discussed in 


Chapter IV. 


II. MODELS AND METHODS 


The model consists of two parts. One is a statistical 
model and the other is a cost model. The type of statistical 
model has been recognized early and studied by many under the 
name of changepoint problem (Carlin, Gelfand and Smith 1992). 
The combination of the changepoint problem with a cost model 
to estimate future costs to make a decision is novel. In this 


chapter we describe the statistical model and the cost model. 


A. STATISTICAL MODEL 

We consider a model for a measure of subsystem performance 
in successive periods of time. Two such measures) of 
performance are the mean number of failures in a time period 
or the mean number of maintenance actions in a time period. 

Consider a sequence of random variables with the following 
structure: 

1. X,, Xz,..., Xo are identically and independently 

distributed, while 

Oe Xe,4, Xc,2,7...- X,wexmabitwa linear trend. 

The time of onset of subsystem degradation, C, called the 
changepoint, will realistically be unknown, as will the 


Magnitude of the linear trend. We will assume 


X, ~ N(p,07), O<isC; 
(2.1) 
~ N(w+(i-C)n,o7), Crisi. 

This is shorthand for the assumption that xX; is 
normally/Gaussianly distributed with mean uw and variance o? up 
to the changepoint time C and is normally distributed 
thereafter, but with mean that grows (with the slope, 1) 
linearly thereafter; wW is the mean number of failures (or 
maintenance actions) in each time period before the onset of 
degradation; C is the time of onset of degradation; TN is the 
slope of the linear trend after degradation; and the variance 
Oo? is a measure of the variability of the actual number of 
failures (or maintenance actions) about mie true mean. This 
model should be appropriate for subsystems whose mean 
failure/maintenance rate per time period, e.g., month, is 
reasonably large, but whose variance is relatively unchanged 
when and if a change in the mean occurs. 

Figure 1 shows data simulated from the above model with 
mean w=4, slope N=1.5, variance o*=1, and changepoint C=10. 
The normal random variables are generated using the Box-Muller 
technique (see G.S. Fishman 1978) which is described in 
Appendix A. 

In Figure 1 the x’s represent the actual data (e.g., mean 
number of failures or mean number of maintenance actions in a 
month), while the dashed lines represent the true, but hidden 


trend. 


NO. of Failures Expected number of Failures 
(Maintenance Act ions) Expected number of Maintenence Act ions 





V-Axis scale 10 see Press <Enter> to exit «ns 
Figure 1. Simulated data (w=4,n=1.5,0%=1,C=10). 

Tt is plain that there is little evidence of any change in 
the demand level until time t=12 at the earliest, when a 
retrospective look suggests that a change took place at time 
t=10. Successively more confirmation is given by observations 
after t=14. If the trend continues as suggested, greater and 
greater confirmation of its direction and magnitude becomes 
available; this can be quantified by the statistical methods 
described in Appendix B (maximum likelihood method) or 
Appendix C (Bayesian method). The statistical methods provide 
estimates of the true trend (denoted by the dashed line above) 
based on the mean number of failures or mean number of 


maintenance actions observed (the x’s above); the estimates 


are of the time of onset of subsystem degradation and the 


magnitude of that degradation. 


B. COST MODEL 

The following is a brief discussion of the cost model 
which will have as input the estimates of the time of onset of 
subsystem degradation and the magnitude of the trend. The cost 
model includes a fixed cost and schedule for upgrading the 
subsystem as well as costs for each failure, maintenance 
action, and AV-DLR action for both the current and upgraded 
subsystems. The detailed cost model which depends on these 
estimates 1s presented in Section B of Appendix D. The model 
is briefly described in the following. 5 

There are three basic unit costs incurred by the current 
subsystem. There is a cost cop (respectively c,,) incurred each 
time the subsystem fails, (respectively requires a maintenance 
action). There is also a cost of qc, for each AV-DLR action. 
These unit costs are computed from cost data in the ROI 
program. The details of their calculation appear in Section A 
of Appendix D. te: 

The proposed upgraded subsystem also has unit costs Cyp 
(respectively cy) incurred each time the subsystem fails, 
(respectively requires a maintenance action). There is also a 


cost of cy, for each AV-DLR action. Details of their 


calculation also appear in Section A of Appendix D. 


An additional cost is an initial fixed cost for the 
upgrade, c,. In addition to the above costs, there is a 
required lead time L to prepare for the upgrade. There is also 
a period of time J to install the upgrade; this period of time 
depends on the installation rate. Details of their calculation 


appear in Section C of Appendix D. 


C. ASSESSING UNCERTAINTY 

Since the estimates of the true trend have variability, 
the estimated future cost of the current subsystem will also 
have variability. As with the estimated future cost of the 
true trend, one can expect the estimates of the future 
subsystem cost to be quite variable until sometime after the 
onset of degradation; this variability is due to uncertainty 
in the estimates of the true trend. 

It is important to consider this uncertainty in the 
assessment of whether or not to upgrade the current subsystem. 
For example, it may be that the estimated mean future cost of 
the current subsystem is larger than that for the upgraded 
subsystem but that the uncertainty associated with the 
estimated mean future cost of the current subsystem is high. 
This may indicate that it is better to wait to accumulate more 
information concerning the apparent degrading trend before 
deciding to invest in the upgrade. Procedures to assess the 
variability of the estimated cost advantage of upgrading the 


subsystem (future mean cost of the current subsystem minus 


future mean cost of the upgraded subsystem) are « 


Appendix E. 
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IIIT. INPUT DATA 


The input data for a radar transmitter on the F-14A that 
are used in the ROI decision aid are attached in Appendix F. 
The data come from the NALDA data base. The model developed in 
this thesis and implemented in the FCT.PAS program is designed 
to use these data to provide a rational time at which to 
upgrade the subsystem, given current information. The main 
difference between the ROI decision aid and the procedures 
developed here is in the forecasting of the future performance 
of the current subsystem. The ROI procedure uses the average 
values of performance measures for the last 24 months to 
forecast the current subsystem’s future behavior. The model in 
this thesis estimates a time of onset of subsystem degradation 
and the magnitude of the degradation using the time series of 
the measures of performance. 

The two time series considered are obtained as follows. 
Table 2 in Appendix F gives the total number of aircraft and 
the total flight hours per month. Because the total number of 
aircraft and total flight hours per month change over time, we 
use the mean number of failures instead of the real number of 
failures in each month to estimate the performance of the 
subsystem. This value is the mean monthly flight hours, 
computed as the average number of systems times the average 


use per month appearing in Table 1, divided by the mean flight 


ee 


hours between failure (column MFHBF in Table 3) for each 
month. The average number of systems and the average use per 
month use the average values of measures for the last 24 
months in Table 2. The mean number of maintenance actions each 
month is estimated as the same mean monthly flight hours as 
above divided by the mean flight hours between maintenance 
action (column MFHBMA in Table 3). 

The number of AV-DLR actions per month depends on the 
number of BCM’s in each month. Examination of the data 
suggests that the number of BCM’s per month is independent of 
the number of flight hours in that month. Thus, we assume the 
number of AV-DLR actions has a constant rate per month. We use 
the average value of BCM’s for the last 24 months to estimate 
the AV-DLR action rate per month. This value is the same as 
that used in the ROI decision aid. The detailed description of 
the parameters used in our decision aid is presented in 


Appendix D. 


eZ 


IV. RESULTS 


Appendix G presents the results of using our decision aid. 
Section A presents results of the procedures using simulated 
data. Section B presents results of using the procedures on 


data for a radar transmitter on the F-14A. 


A. RESULTS FOR SIMULATED DATA 

We first discuss the results for a set of simulated data. 
Two time series using data simulated from model (2.1) were 
used. Data for the mean number of failures in each month was 
simulated from the model with parameters u,=225, N,p=4, 6,’=144, 
and C,=25. Data for the mean number of maintenance actions was 
Simulated from the model with parameters u,=440, n,=5, 6,’=225, 
and Cy=25. The length of both time series is chosen to be 40 
months and the time horizon is to be 150 months. These data 
are chosen to approximately mimic the real data discussed in 
the next section, but to have more apparent linear trends. The 
other parameters concerning cost computation were chosen to be 
equal to those in the next section which also gives a detailed 
explanation of the computation of the parameters. Figure 2 in 
Appendix G displays graphs of the two sets of simulated data 
along with the true mean as a dashed line. A listing of the 


Simulated data appears below Figure 2. 


3 


For each time 6St<40 the maximum likelihood procedure and 
Bayesian procedure were used to estimate the model parameters 
using data x,)9)x,/ eee 

The cost of upgrading at time tein theefuture C,(epe) ane 
the cost of never upgrading Cagicy) were computed using the 
expressions (D.8) through (D.13) in Appendix D. 

The parameters used in the cost model appear in the menu 
below the listing of the simulated data in Appendix G. A 
description of how they are used appears in Section A of 
Appendix D. The time t which minimizes €,(t,t) is computed. If 
C(t) < min, C,(t,t¢), then the minimizing tT as taken to wee ame 
horizon time, taken to be 150 in the example. A positive value 
of the cost advantage of upgrade, E(t) -Gy(t, t), indicates that 
it is better to switch to the new subsystem; the more 
positive, the greater the estimated advantage of changing to 
the new subsystem. 

Below the listing of the menu in subsection 2 of Appendix 
Gis a listing of the minimizing tT for each time t for the 
Bayesian procedure. Also displayed is the mean cost advantage 
of upgrade for the best policy and its two standard deviation 
bounds for assessing uncertainty; if the minimum cost policy 
is never to upgrade, then the cost advantage is C,(t)-C,y(0,t). 
Notice that initially the best policy is to do nothing. 
However, as time increases, the best policy is to begin the 
procedure to upgrade immediately. This policy is interspersed 


with the policy to do nothing until time 26. After time 26, 


14 


the best policy is always to begin upgrading immediately. The 
true best policy would be to do nothing up to time 26 and then 
to start upgrading immediately. Thus the policy using the 
estimated costs occasionally gives false alarms, (suggesting 
that the upgrading start before time 26). This behavior 
Suggests that it is prudent to wait for confirmation of a 
decision to upgrade before starting the upgrade policy. 
Figure 3 presents graphs of the cost advantage of 
upgrading obtained from the Bayesian assessment of cost 
VMamtability for minimum cost policies from t=6 to t=40 as 
described in Section B of Appendix E. The graph displays the 
mean cost advantage of the minimum cost policy and the mean 
cost advantage plus and minus two standard deviations; if the 
minimum cost policy iS never to upgrade, then the cost 
advantage is C,(t)-C,(0,t). The width between two bounds can 
be interpreted as representing an approximate Bayesian 
posterior density for the true expected or mean cost 
advantage, given observations up to time t. The width between 
the bounds becomes smaller as more data accumulates and the 
uncertainty of the estimates of the changepoint and the 
degradation rate, 1, is reduced. The width of the bounds 
provide prospective on the risk of changing soon, or waiting. 
Apparently the chance of making the wrong decision decreases 
1f the decision maker waits, but also the value of making the 
more nearly correct decision decreases, for there is less time 


to the horizon. Recall that the true time of onset of 


5 


subsystem degradation occurs at time 25. If one waited until 
the lower confidence bound of the cost advantage of the 
minimum cost policy becomes positive, then one would wait 
until time 34 to make a decision to upgrade the subsystem. 

Subsection 3 presents the output for the likelihood 
procedure using the same simulated data. The best policy is 
always to upgrade immediately after time 26. Comparing the 
upgrading policies with those from Bayesian procedure, the 
only difference occurs at time 25; the Bayesian procedure does 
not recommend upgrading, while the likelihood procedure 
recommends upgrading immediately. This is before the time of 
onset of degradation. Figure 4 presents graphs obtained from 
a bootstrap assessment of the variability of the cost 
advantage of the best policies from time 6 to time 60. The 
variability of the cost advantage of upgrading estimated from 
the aaNENIeTE procedure tends to be larger than that from the 
Bayesian procedure. This result is due to the likelihood 
procedure’s having more variability in the estimated time of 
onset of the degradation. 

Simulation was also used to investigate the behavior of 
the estimation procedures for two other changepoint models. In 


one (a jump model) 


16 


X, ~ N(p,0*) ore, 
(4.1) 
~ N(p+,07) C+isi; 
that is, X; 1s normally/Gaussianly distributed with mean p and 
variance O* up to the changepoint time C and is normally 
distributed thereafter with an different constant mean +6. 
The other model (jump plus linear trend model) 
X, ~ N(p,07) lanier 
(4.2) 
~ N(p+d+n(1-C),07) C+1si; 
that is, X; is normally/Gaussianly distributed with mean uw and 
variance oO” up to the changepoint time C and is normally 
Gaoerwbutedsthereafter with a mean that has a jump 6 and then 
grows (with the slope N) linearly thereafter. 

Simulation studies indicate that the maximum likelihood 
procedures for these last two models yield estimates that are 
more sensitive to the local behavior of the data than those 
for the model presented in Chapter II. This sensitivity to 
local behavior tends to produce more "false alarms" concerning 
thempresence GfWdegradativonweThevestimator» for 6 in both of 
the above models is also greatly influenced by local behavior 


in the data leading to a large assessment of variability. 


B. RESULTS FOR DATA FROM A RADAR TRANSMITTER ON THE F-14A 
Results concerning the analysis of data from a radar 


transmitter on the F-14A appear in Section B of Appendix G. 
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The mean numbers of failures each month, and the mean 
numbers of maintenance actions each month, were calculated as 
described in Chapter III. A listing of the two time series 
appears below Figure 5 which presents graphs of the two 
series. The data show considerable variation and tend to 
increase. Thus the subsystem appears to be degrading over 
time. There is no way of knowing what actions, if any, were 
taken in the event that such a tendency was noted. Using the 
parameters calculated in Appendix F, the cost model for the 
example has the following features. The fixed cost per failure 
for the current subsystem is cC,,=1,323.6. The fixed cost per 
maintenance action iS C,,=166.8; and the fixed cost per AV-DLR 
action iS C,,=1,120.3. The mean number of AV-DLR actions per 
month is 7.8. There is a fixed initial cost c,=5,400,000 for 
upgrading the subsystem. All costs are in dollars. The 
upgraded subsystem is assumed to have a mean number of 
failures of 125 per month, and a mean number of maintenance 
actions of 485 per month. The fixed cost per failure for the 
upgraded subsystem is assumed to be cy-=1,323.6. The fixed cost 
per maintenance action for the upgraded subsystem 1S Cy=204.4; 
and the fixed cost per AV-DLR action for the upgraded 
subsystem, Cy, 1S the same as current subsystem. The 
procedures used to compute the cost parameters appearing in 
this section are described in Appendices D and F. The data 
were collected over 5 years (January 1987 through December 


1991). There is a time horizon, H=180 months, during which 


rs 


this subsystem or its upgrade will be used, and a lead time, 
L=60 months, to prepare the upgrade. The decision to upgrade 
depends on the estimates of the time of onset of subsystem 
degradation and of the magnitude of the trend. The assessment 
of the cost of upgrading should reflect the uncertainty of 
these estimates. 

For each time t26, the following policies are considered: 
upgrade the subsystem at each future time until the time 
horizon; all potential upgrading times from the present time 
until the time horizon H-L are considered; that is, if the 
current time is t=40 then the policies that would upgrade the 
subsystem at time 40, time 41,..., time 120, (which is H-L), 
are considered. For each current time t, the (estimated) costs 
of these policies are compared to the (estimated) costs of 
never upgrading the subsystem. The "optimal" (maximum 
estimated cost advantage) policy can then be found. 

For each time ¢t26, the model considers the data 
accumulated up to time t and using the data as of that time 
estimates the time of onset of subsystem degradation and 
magnitude of the trend. The time of onset and the magnitude of 
the trend are estimated for each of the two time series 
(numbers of failures and numbers of maintenance actions) 
independently. For each current time t, the estimated mean 
cost for each policy to upgrade the subsystem at some future 


time iS computed using the current estimates of the trend. 
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The column "Best upgrade time" in Appendix G presents the 
times to upgrade the subsystem which correspond to the maximum 
estimated mean cost advantage policies for each current time. 
If the maximum estimated mean cost advantage policy is never 
to upgrade (negative mean cost advantage), then the time to 
upgrade is set equal to the horizon time, H=180. These results 
are shown on the computer screen or contained in the output 
file. Also displayed is the estimated mean cost advantage of 
the best policy, standard deviation of the cost advantage, and 
the mean plus or minus two standard deviation bound for the 
best policy for each time t. If the best policy is never to 
upgrade, then the mean and standard deviation of the cost 
advantage are computed for the policy that starts to upgrade 
immediately. 

We first describe the results obtained by using the 
Bayesian procedure. These results appear in subsection 2 of 
Section B. The best policy for current times 5-20 are either 
to upgrade immediately or to never upgrade. We can compare the 
data with the policies. The data vary somewhat before time 20, 
but there is no evidence of a trend occurring. Because of the 
variation of the data, the results are unstable. An apparent 
trend appears after time 20. The Bayesian procedure suggests 
this subsystem should be upgraded immediately for all the 
times after time 20 except time 46 and time 50. The declining 
data values at times 40 to 53 appear not to be able to 


overcome the increasing numbers before them. We can also check 
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the "cost advantage of upgrade“ column. Even these costs can 
not accurately represent the real future cost advantage, but 
it still gives further information for the subsystem upgrade. 
At time 46 and time 50 the small negative mean cost advantage 
compared with others do not give strong evidence to maintain 
current subsystem. Figure 6 presents graphs obtained from the 
Bayesian assessment of the variability of the estimated mean 
cost advantage for the best policies from times 6 to 60 as 
described in Section B of Appendix E. 

The results of the maximum likelihood procedure are 
presented in subsection 3. We can compare these results with 
those of the Bayesian procedure. The only difference in best 
policies occurs at time 50. It changes from never upgrade to 
upgrade immediately. We also compare the estimated mean cost 
advantage of the best policies with those computed by the 
Bayesian procedure. Again, the variability of the estimated 
mean cost advantage is larger than the variability in Bayesian 
procedure. Figure 7 presents graphs obtained from the 
bootstrap assessment of variability obtained by the maximum 
likelihood procedure of the cost advantage of upgrade ¢,(t)- 
Gaye) for the best policies from t=6 to t=60; if the best 
policy at time t is never to upgrade, the cost advantage of 
the best policy is €,(t)-C,(0,t). The numerical values of mean 
cost advantage and two standard deviation bounds appear after 


the best policies. Displayed is the mean cost advantage 


Zr 


100 
- = 1 A s -f 8 
m(B;%),60) = aioajpeaeoue: o2) Cy(b5t,60)], (4.3) 


and the mean plus and minus two standard deviations, where T, 
is that time which maximizes the cost advantage for the 
original data; if the best policy is never to upgrade, then 


To=0. The bootstrap variance is 


100 


1 n a 
E2(B;t,,60) = 99 pz [0 PUG e Degree taeeteren ten (4.4) 


Displayed is 


m(B3t,,60)e2 ce (ar u., OO (4.5) 


Comparison of Figures 6 and 7 indicates that the bootstrap 
estimates of the variability of the estimated mean cost 
advantage of the best policies is larger than those for the 
Bayesian procedure. This larger variability for the maximum 
likelihood procedure is due to more variability in the 
bootstrap distribution of the estimated time of onset of 
subsystem degradation. The maximum likelihood procedure 
appears to be more sensitive to local features in the data 
than the Bayesian procedure. 

The two standard deviation bounds appearing in Figure 6 
and Figure 7 indicate that the variability of the estimated 
mean cost advantage becomes relatively small and stable after 
time 30. Notice that the estimated mean cost advantage is 


positive, which suggests that it 1s advantageous to start an 


Za 


upgrading process. However, since the lower confidence bound 
is negative, it may still be worthwhile to wait for more 


evidence before starting the upgrading program. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Preliminary mathematical models have been formulated for 
the possible onset and growth of subsystem degradation. The 
model recognizes that the time of onset of a degrading trend 
may be random, and hence initially unknown, and that the trend 
Magnitude is also initially unknown. The trend magnitude will 
become better known as more data is accumulated. Statistical 
procedures have been developed to estimate the time of onset 
and the trend magnitude. A cost model that is compatible with 
the existing decision aid, the ROI procedure, has been used to 
develop procedures (which recognize the uncertainty concerning 
the time of onset and magnitude) to determine estimated costs 
and the associated risks of upgrading the subsystem at 
different times in the future. An experiment uSing real data 
gives reasonable results and indicates that the consideration 
of variability in policy costs due to uncertainty concerning 
the time of onset and trend magnitude can lead to wiser 


decisions. 


B. RECOMMENDATIONS 
Two procedures to estimate the mean cost advantage and 
assess its variability are considered. One uses maximum 


likelihood for estimation and the bootstrap to assess 
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variability. The other uses a Bayesian model. The Bayesian 
procedure requires much less computational effort than the 
maximum likelihood/bootstrap procedure and appears to give 
Similar results. Thus we suggest that the Bayesian procedure 
be used to estimate best time to upgrade and to assess 
variability of a cost advantage. 

The changepoint model considered in this thesis has a 
linear trend after the changepoint. A linear trend may 
overestimate the magnitude of the degradation. Other 
possibilities exist. For example, another possible model is 
that the trend be proportional to the square root or some 
other power less than 1 of the time since the changepoint; 
Ehat is, 

X, ~ N(p,o7), O<i<C; 
(Ste) 
~ N(pt/i-Cn,o7), Crisi. 
Future work can extend the estimation procedures to such 
cases, and study the sensitivity of change policies and their 
costs to different specifications of degradation growth. 

The cost of the subsystem in this thesis is the sum of 
costs due to failures and maintenance actions. In this thesis 
the costs due to failures and the costs due to maintenance 
actions are estimated separately. Future work can extend the 
estimation procedures to multivariate time series, if 


appropriate. 
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APPENDIX A 


NORMAL RANDOM NUMBER GENERATION 


A simple scheme for generating Normal random variables is 
the Box-Muller technique (G.S. Fishman 1978). The procedure 
generates two independent standardized Normally distributed 
variables X and Y as follows. 

1. Generate U,, U, as independent random variables 

uniformly distrabutedtone 07a 


2. Set 


X = /-2(InuU,) cos (2nU,)- 
(A.1) 


Y = /-20@nv|) sim.) 
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APPENDIX B 


MAXIMUM LIKELIHOOD ESTIMATION 


Suppose observations of the variables (numbers of failures 
or maintenance actions during time periods 1,...,t) X;, X5,..., 
X, are available; denote them by x,, X,..., xX,. Then the 
likelihood function for the unknown parameters, I, C, TN, 0’ is 
as follows for the model in Chapter II (Donald P. Gaver and 
Patricia A. Jacobs 1992). Since the number of failures 
(maintenance actions) in successive time periods are assumed 
to be independent, for time 1 <CS# t¢t, the likelihood function 


is 


Cc - (xy-p) 2/20? t - (x ,-p- (1-C) n) 2/202 
L(p, C, n, 02; data) jee —_—_— Il ee (B.1) 
1=1 v2no" 4=C+1 ¥2no" 
so the log-likelihood is 
Cc 2 t 4 2 
Sas -(x,-p-(i-c 
1(p,C,n,07; data) aighpy ia + 2 ee 
1=1 207 1=C+l1 20 (B 2) 
- = Ino? + constant. 
This can be concisely written as 
t ms ee i-c)’ 2 
l(p,C,n,07; data) = > aa 
“9 (B.3) 


- =1no? + constant. 


Ze. 


where 


1-C ifeiI>oC, and 
(7=CcCy = (B.4) 
0 LE ode 


Note that the above applies if there is a changepoint within 


the range of observation; otherwise, if C > t then 


c ~ (x4-p) 2/20? 


L(p,C,n, 07; data)  ————— (B.5) 
1=1 270 
and 
: =a) t 
l(p,C,n,07;data) =} ee 5 ino? + constant. (B.6) 
1=1 0 


Now in the following hold C fixed and behave as if it were 
known and the objective is to maximize 1 with respect tou, 1, 


and 0”. Begin by differentiating with respect to uw: 








=n : sie (1-6 
Ol - Se aS Xe eee 
ow fi (o? 4=C+1 o? 
(B.7) 
t 
=> “iP if ot. 
1 «(o* 
These expressions can be simplified and combined: 
E=C 
o2 91 . ex(t) = ti —ippee 
be jJ=1 
(t-C)2+ (t-C) (B.8) 


tx) — "ei -{ 5 In. fOr t>C 


tx(t) -“Gaye for t<c 
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where 


(B.9) 
i=1 
Rewrite this as 
IR oe (B.10) 
a t(x(t) -p) 5 n 
where 
t-C if t2c; 
ai (3.11) 
0 TF sac. 


If the derivative is set equal to zero we obtain the first 


"normal equation" 


Bost y, (Cc, t)n = x(t) 


(B.12) 
where here 
Seale. Cet. (C=C) (B.13) 
wp, (Cc, €) | 5 |, 


Next differentiate (B.3) with respect to Tj: 


92.92 = (x,-p-(4-Cy'n) (4-C 
om 4 H n 


2 =7 


zc t t (B.14) 
pee yy (2c) — 9 y ((1-c)")? 
a7 i=1 


GX5(G, Gey Gun (C, C) + nity, (Cc, t) 
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where 


XG.e) = 2) x,(4- ae (B.15) 
=1 


W2(C, t) = 4 (d= cy)? = zeae (he ON Leet oaea (B.16) 
CTH 


Set the derivative equal to zero to obtain the second normal 


equation 


Y, (Cre 9b (Chey n = XC, ec). (B.17) 


Dreferentiate wien respect scouG. 


Ol 


C 
1( ‘ a Ee B.18 
= = Diy islet ye =)? 5 ta) ( ) 


if this is set equal to zero and solved for o° there resi 
1 C 
68..= Ee (x;=8 al 4C)h))* (B- 


Now solve the first two normal equations for the maximum 


likelihood estimate, conditional on Cc; the resmieers. 


_ Box = x 

fi(c) = +. - ee = ty)? (B.20) 
_ = - Wx 

Rica my? y.)2 (B.21) 


for C<t; for C2t, fi(C)=x, f\(C)=0. These can now be substituted 
into (B.19) to obtain the maximum likelihood estimate for 0’ 


in terms of the other estimates, all conditional on the value 
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of Cc. Finally substitute the above estimates into the 


expression for the negative of the log likelihood: 


S(C;data) = -21(A(c, t),C,(C, t) ,82(C, t) ; data) 


2 lpm (a7hl(Ct))? 9 1 GQ (xi-A(C, t) -(i-c) (Cc, 0)? 
eh; ‘sz: mae 


82 (C, t) 0-(C,E) 
+ 1n@2(C, t) 


1. + 1n62{ Ct) (B.22) 


and obtain the value of C that minimizes S(C;data) over the 
lege (il, 2,..., t€); denote this by C(t); the last equality in 
the above expression follows from the definition of 6°(C,t) 
given by (B.19). Thus, the estimate of C is chosen to minimize 
the sum of the squared residuals. If the minimum of S(C;data) 
occurs at t=C, then the conclusion is that no change has 
occurred in [1,t]. Note that all estimated parameter values, 
namely , fl, and 6* depend upon the C value in use, and so the 
dependence of S upon C involves that implicit dependency. Once 
C(t) is developed this value is substituted into the 
expressions for fi, fj, and 6’ to obtain the maximum likelihood 


estimates of those parameters. 
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APPENDIX C 


BAYESIAN ESTIMATION 


An enhanced version of the basic model presented in 
Appendix B is obtained by assuming that the changepoint (time 
of onset of degradation) is a random variable, C, with 
specified distribution whose parameter is unknown and subject 


to a probability density, "(-). Specifically, suppose 


Pic=kj) = (1-p)‘"p, (c.1) 
1.e€. 1S geometric, and that the parameter p has a prior 
distribution, Tie 

We also use the linear normal model here. Putting 
(uninformative) priors on u, N, and p, it is shown that the 
joint posterior density of those is’ straightforwardly 
obtained; the parameter oO? is initially estimated from 
residuals. In principle all of the above could be carried for 
any arbitrary, but reasonable, discrete distribution that 
might better represent what is known about the changepoint 
process. A similar statistical model was used by Smith (1975). 

In what follows we sketch the development (Donald P. Gaver 
and Patricia A. Jacobs 1992). Suppose that observations, 


X1,..., &, are available up to time t, it follows that 


S2 


P{pe (dp) , C=k, pe (dp) , ne (dn) |X,=x,, ...,Xy=%y,.. . » X,=X_, 97} 





c 
n(p)(1-p)k"*p [J Z explo (x,-#- -kI'n)?}dudndp for kst 


i=1 §/2n02 20° 








t 
n(p) (1-p)*"2 z exp 
I] ¥27o7 


ee (-3¢-1)?}dudn dp for k=t+1 (C.2) 


where 


ls 


i-k if izk, and 
(c.3) 


0 if 1<k. 
The term involving (1-p)‘ represents the case in which no 
changepoint has occurred; we will set k=t+1 for this case. 
By a completion-of-squares process one can write the 
likelihood function for given C=k as a bivariate normal 
density with parameters dependent on k and data up to t; the 


exponential term of the likelihood is written as 


{joo 


4=1 


a (aw kya) | (C.4) 





= cm |e =H)" -2p (Heb) inn) seuoen) . (n=) © a?) - tk, of 
2 (1-p?) Y? YV v2 


for 1S$kst-1; for k2t we have no changepoint so the exponential 


term of: the likelihood is of the form 


Bi3 





(c.5) 
=o exp 5 eB -K(t, e| 


where in the above c is a constant, and the parameters all 
depend upon k, t, and x(t), the data up to time t. 
For k<t, the parameters of the bivariate normal (C.4) turn 


out to be 


Xe) Y, oa Xo (k, E) WY, 








p(k, t) = 7 an: : (C.6) 
i(k, t) = a0 eye (C.7) 
y?(k, t) = a ee S (c.8) 
v2(k, t) = eee (c.9) 
and 
p(k,t) = - mae (C.10) 
Vv. 
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where 


— all 2 
x(t) ==)o xy; 
Cia 


c 
aif mt si ate 
1 c 
y, = = isk); 


1=1 


es 
v, ey): 


and 


c 


ae tee < Y> (x,-B(k, t) “Wk, t)(i-k)); 


1=1 


For the case ket 





and 





1=1 


BS 


(C.11) 


(C.12) 


(C.13) 


(C.14) 


(C.15) 


(C.16) 


(C.17) 


(c.18) 


N(k,t)=0, v*(k,t)=0, and p(k,t)=0. These values can be derived 
directly from (C.4) and (C.5); the procedure is similar to 
that in Appendix B. 

If the bivariate normal form is utilized in (C.2) and the 
integration is performed over p we obtain the joint 
conditional density of C, U, and N given the data and o’ in 


the form 


P{C=k, pe (dp) ,n€ (dn) |x(t) , 07} = 


n*(k, t) en 1 (L2H)? 2p Leena), (n-A)*)|(¢,19) 
2nV1-p?yv 2(1-p*)\_ oy’ in 


for k<t where 


m*(k,t) = ele (1-p)*-? px(p) dpexp(-K(k, t))(2nV1-p?yv);  (c.20) 





fom kee 
=k, ,02\ = w*(k, t) — -1_ (u-p)? C.21 
P{c=k, pe (dp) |x(t) , 07} = n*( ) _ere| ain Ht) ( ) 
with 
n*(t,t) =c* im (1-p)*" px(p) dpexp(-K(t, t)) (V2ty) (Cc. 22) 
m* (21,5) — otf (1-p)'x(p) dpexp(-K(t+1,t))(v2ry) (Cc. 23) 
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and 


eal ae 
» I, (1-p)*"" px(p) dpexp(-K(k, t)) (2ny1-p?yv) 


k=1 


+ fo a-pi* pap) dpexp(-K(t, t)) (Vary) (C.24) 


+ [0 (1-p)*n(p) dp exp(-K( +1, €)) (VaR) 


Note that {x (k,t), kst} is the marginal probability that the 
changepoint occurs at any time k up to and including t; while 
m (t+1,t) is the posterior probability that no changepoint has 
occurred up to time t. 

For each time t, the estimate of o” is computed from the 
squared residuals for each possible value of C=k in the 


following manner; let 


t 


= 1 — = Dean : , 
Gc, ©) = ieee ig Xe Back eiok) | 2eeRst; (C.25) 
c 
Cee ee weket) |? YF koe. (¢.26) 
t-1im 


Finally, the estimate of the variance o* based on data x,,... 


x, ass 
t+1 
67(t) = d) 1° (k, SO -(k, 6) (€.27) 
k=1 
Given C=k, k<t, and the data x,,..., X,, the posterior 
dastribution of (u,T)) is bivariate normal with mean 


(u(k,t),N(k,t)), variance of w equal to y?(k,t), variance of 


oF] 


N equal to v*(k,t), and correlation p(k,t); for k=t, t+1, =0 
and the posterior distribution of w is normal with mean U(k,t) 
and variance y’(k,t). Hence, given the data x,,..., x,, the 
posterior distribution of (uw,N) is a mixture of bivariate 
normal distributions with mixture distribution {nx'(k,t), 
kSt+1}. 

Since the bivariate normal has 5 parameters to be 
estimated, the estimation procedure begins with data x,,..., 


x,. The initial estimate of o? is 
3 1 = C.28 
° _s\2 ; 
7 ae (x, x) ( ) 


where x is the sample average of the first 5 data points. For 
each time t, estimates of the posterior distribution are 
obtained from equations (C.6)-(C.26). The updated estimate of 
Oo? is used as input for the calculations for the next time 


period. 
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APPENDIX D 


A MATHEMATICAL MODEL FOR COST AND RETURN-ON-INVESTMENT 


In this Appendix we describe the cost model. The cost 
model is formulated to reflect the costs that are used in the 
decision aid in current use, the ROI procedure. 

There are several different unit costs. There is a cost 
due to subsystem failures; a cost due to subsystem maintenance 
@e=ion; and a cost adue"to AV=DLRwaction”. 

In addition there is a planned horizon H during which the 
parent system will be operative; when the horizon is reached 
all (remaining) parents are stored or disposed of. There is 
also a lead time L before the upgrade is initiated and an 
installation period of length J. 

Section A below describes the calculation of the unit 
costs. Section B below describes the cost model in detail. 
Section C describes the cost estimation procedure. 

Let C,(t,t) denote the estimated mean cost of deciding at 
time t to begin the procedure to upgrade the subsystem T time 
Ungitcman thewtutwune. Let €,(tJedenote the cost of deciding 


never to upgrade. If 


min Gobet) < Celt) (D.1) 


then it may be advantageous to begin upgrading process at that 


time tT which minimizes the left hand side of (D.1). However, 


eg 


a large amount of uncertainty concerning C,(t,t) and C,(t) may 
indicate that it 1s better to postpone the decision to upgrade 
until more data has been obtained. Appendix E describes 


procedures to assess the uncertainty of the cost estimates. 


A. UNIT COSTS 

This section describes the computation of the unit costs. 
All the parameters used in the computation of cost can be 
obtained from the ROI program. These costs for a radar 
transmitter for the F-14A can be found in the Appendix F. 
There are 5 costs in the ROI procedure for both the current 
and upgraded subsystem. The 5 costs are computed as follows 

1. Using the "O" MH/MA times the MA at "O" level (MA/ML1 
Ratio) gives the total number of hours ("O" MH) spent by the 
squadrons on maintaining this subsystem. Multiplying the "O" 
MH by the composite rate for "O" level (ML1 Rate) gives the 
“O" level manpower cost. 

2. Similarly, using the "I" MH/MA times the MA at "I" 
level (MA/ML2 Ratio) gives the total number of hours ("I" MH) 
spent by the squadrons on maintaining this subsystem. 
Multiplying the "I" MH by the composite rate for "I" level 
(ML2 Rate) gives the "I" level manpower cost. 

3. Again using the simple relationship of MH/F times the 
VF gives the total number of hours spend on repairing this 
subsystem. Multiplying the hours by the composite rate for "I" 


level (ML2 Rate) gives the manpower cost per repair. 
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4. The cost of materials used for repairs is the product 
pEmCoseMOreinaterialie™per repair (S/YRpr) times the number of 
repairs. 

5. The total AVDLR cost is the product of the number of 
BCM systems times the AVDLR unit cost. 

The calculation of the model unit cost parameters due to 
failure, maintenance action, and AV-DLR action is summarized 


as follows 


Cow = (“O" MH/MA) 01¢(MA/ML1 Ratio) (ML1 Rate) (D.2) 
+ ("I" MH/MA) 514 (MA/ML2 Ratio) (ML2 Rate), : 

Cu = (“O" MH/MA) 1¢,(MA/ML1 Ratio) (ML1 Rate) (D.3) 
+ ("I" MH/MA) joy (MA/ML2 Ratio) (ML2 Rate), 

Cor = (MH/F) o1q(ML2 Rate) + ($/RDr) oig, (D.4) 

Cur = (MH/F) poy (ML2 Rate) + ($/RPI) ney, (D.5) 

Co, = (AVDLR Cost/Unit) oj, (D.6) 

Cua = (AVDLR Cost/Unit) ney: (D.7) 


The subscript O represents current subsystem; N represents 
upgraded subsystem; F represents failure; M represents 


maintenance action; and A represents AV-DLR action. 


B. ESTIMATED FUTURE MEAN COST 
Fix a time t and let ¢, (respectively C,) be the estimate 


obtained at that time of the time of onset of subsystem 
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degradation due to failures (respectively maintenance 
actions); let fl, (respectively fl,) be the estimate of the 
constant mean number of failures (respectively maintenance 
actions) in each time period before the onset of subsystem 
degradation; and let fj, (respectively f\,) be the estimate of 
the magnitude of the linear degrading trend in the mean number 
of failures (respectively maintenance actions) after the onset 
of subsystem degradation. 

The estimated future mean (total, undiscounted) cost due 
to subsystem failures of a policy that initiates upgrading T 
time units in the future is 


J 
Corll » | (+ L41) +>) U1 -ans) J 


s=1 





+ CypAp 


S=1 





J 
Fas) tH (ersten) if Coe 


CA cree) — 


» (fi ptf (st (t-E,)*) ) 
s=0 
Cor Ltt+T 
+ SS (@ptfip(st(t-C,)*) )(1-a (s-(t+L))) 
S=Ltt+1 
J 
+ Cype | >) @(S) 1 H—( Cer) LLCESG (D.8) 
s=1 





where L is the lead time to begin installation of the upgraded 
subsystem; A, is the mean number of failures per month for the 
upgraded subsystem; J is the length of the installation 


period; a(s) is the fraction of subsystems that have been 
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upgraded s time units into the installation period; and H is 

the time horizon, (the useful lifetime of the subsystem). 
The estimated future mean cost due to. subsystem 

maintenance actions of a policy that switches to the upgraded 


subsystem T time units in the future is 


J 
(t+L+1) +} (1-a(s) ) 


s=1 


Conf 1 








S=1 


J 
“ cach a5) +H (ers+Z-70) ate ate 


~ = L+t pe 
Cig (Ts t) > (ft y(st+(t-C,)*) ) 
s=0 
Com Letts 
+ a (fiutfiy (st (t-€,)°) ) (1-a@ (s-(t+L))) 
SBLt+t+1 
J 
em © 0s) + (H-(trt+L +) ) DCR Es (D.9) 
sl 








The total estimated future mean cost of the policy that 


begins the upgrading process T time units in the future is 


Cymer —uGe tse) = C(t, t) + Cy 


s=1 


J 
+ CaaYa ceszen) +38 (1-a (s) ] (D.10) 
J 
teal yest H-( tec+i4J),) 
s=1 


where c, is the initial fixed cost for the upgrade; yY, is the 
mean number of AV-DLR action per month for current subsystem 


and is computed as the average number of BCM’s for the last 24 
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month; and A, is the assumed mean number of AV-DLR actions per 


month for upgraded subsystem. 


The estimated mean cost due to failures of a policy that 


never upgrades the subsystem is 


Corll p(H- (t-1)) tt Gee 


Conte) (D.11) 


H-t 
Cor |S (Bp+fip(s+(t-G,)"))| if pst. 


g=0 


The estimated cost due to maintenance actions of a policy 


that never upgrades the subsystem is 


Gaui (H=(E-1)) if (ES 
Ca, ( be Per : (D.12) 
Com | >, (Bitty (ee Ce Guwen |) dhteraten 
3=0 


The total estimated mean cost incurred by a policy that 


never upgrades the subsystem is 
&4(t) = Cop(t) + Ca, (6) eereay ae ate) (D.13) 


C. THE ESTIMATION PROCEDURE FOR COSTS 

Table 1 in the Appendix F lists the menu of the ROI 
procedure. All of the cost parameters for the current 
subsystem and all parameters for the upgraded subsystem used 


in our decision aid appear in, or are computed from, the 


numbers in Table 1. 
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The current decision aid, the ROI program, uses the MTBF 
and the MTBMA listed in Table 1; these are obtained by 
averaging the measures of performance for the last 24 months 
in the MFHBF and MFHBMA time-series data which appear in Table 
3. The decision aid developed in this thesis uses all MFHBF 
and MFHBMA time-series data. Our decision aid uses the 
estimated measures of performance for the upgrade subsystem 
appearing in the improved column of Table 1. The cost for the 
upgrade and the length of the lead time are also taken from 
Table 1. The other information section of Table 1 contains the 
cross-over month computed by the ROI decision aid. The cross 
over month is the ROI program’s measure of evaluating the cost 
effectiveness of the upgrade; it is month in which the total 
cost of subsystem upgrade becomes smaller than the estimated 
cost of not upgrading the current subsystem. The other 
information section also contains the number of systems and 
their use per month; these are computed as average values of 
measures for the last 24 months in Table 2 which gives the 
total number of aircraft and the total flight hours per month. 
The other four values in the section are also computed by ROI 
procedure. 

Because the total number of aircraft and total flight 
hours per month change over time, we use the mean number of 
failures instead of the actual number of failures in each 
month to estimate the performance of the subsystem. This value 


1s the mean monthly flight hours, computed as the average 
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number of systems times the average use per month divided by 
the mean flight hours between failures (column MFHBF in Table 
3) for each month. The mean number of maintenance actions each 
month is estimated as the same mean monthly flight hours as 
above divided by the mean flight hours between maintenance 
actions (column MFHBMA in Table 3). 

For each time t the procedure of Appendices B and C are 
used to independently estimate (f,, fY%;, C-) and (f,, fGen 
The parameter Y, 1s obtained by taking average of BCM’s for 
the last 24 months. The number of time periods used for the 
installation, J, 1s obtained as the number of systems in Table 
1 divided by the installation rate (Kits/Mth installed) in the 
same table and rounded up to the next iaieager. This integer is 
then compared to the time remaining (the starting time for the 
upgrade until the time horizon, H-t-L-T), the smaller of these 
two numbers is then used to represent J. The fraction of old 
subsystems that have been upgraded s time periods into the 
installation period, a(s), is obtained as s times the 
installation rate divided by the number of systems. The fixed 
cost to upgrade c, is obtained by adding the cost per kit and 
the cost to install one kit together multiplied by the number 
of systems, and added to the other terms in the section "Cost 
for fix". The lead time L appears in the lead time section. It 
is obtained by adding all the items in the section except 


WRTts/Mthn installea 
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APPENDIX E 


ASSESSMENT OF UNCERTAINTY 


In this Appendix we discuss procedures to assess the 


variability of the cost estimates described in chapter II. 


A. THE BOOTSTRAP PROCEDURE 

A re-sampling technique called the bootstrap can be used 
to assess the variability of the maximum likelihood estimated 
mean cost associated with a policy (B. Efron and R. Tibshirani 
1936) . 

Seon emancd let f(t) ,(t), 6,(t), and C,(t), 
je {F,M} denote the maximum likelihood estimates obtained from 
the data. A bootstrap replication has the following steps. 

1. Using model (2.1) with parameter values equal to the 
Bewamates (i-(t), f(t), 6,(t), and C.(t) simulate data x,,(b), 
Xp (b),..., Xp, (b). Using the simulated data use the maximum 
likelihood procedure to obtain bootstrap estimates fl,(b,t), 
ney et), O-(b,t), amd C.(b,t); b denotes the iP boost rap 
Simulation; b=1, 2,..., B, where B is the number of bootstrap 
samples utilized. 

2. Repeat step 1 for the estimates f,(t), fTly(t), 6,(t), 
endec (eye to sobtain bootstrap estimates f,(b,t), MNy(b,t), 


Ouse), and C,(b,t) . 
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3. Using the bootstrap estimates obtained in steps 1 and 
2 above compute the future mean cost of a policy that switches 
to a new subsystem Tt time units into the future, C,(b;t,t) 
using (D.10). Also compute the cost of never changing C. (bee 
froma(Delaie 


4. Compute the cost advantage of upgrade 


Ci(b;t) - C. (bree (E.1) 
The results reported in Appendix G use B=100 bootstrap 
replications. After the 100 replication are generated, the 
mean and variance of the bootstrap cost advantage are 
computed; that is 


é 


100 
= A tb: t)-G.(b: 
m(B;t,t) = Too py [Coli pcre. oN: (E.2) 


100 
ee A (b:t)-G.(b; -m(B; 
°(B;t,t) = 99 py [oli t)-C,(b;+,t)-m(B;t,t)). “Clee 


B. THE BAYESIAN ASSESSMENT OF UNCERTAINTY 

For series data x,,, Xj2,..., Xj, the Bayesian procedure 
described in Appendix C yields a posterior distribution for 
the time of onset of subsystem degradation as of time t, 
namely C,(t); it also gives estimates of the conditional 
variance of ut,(k,t), Ny(k,t), and their conditional covariange 


given C;(t)=k, for je{F,M}. These estimates together with the 
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cost model can be used to obtain the variance of the cost 
advantage for a policy. 


For example for k<t 


ete ,s(t)-C,.(t, t) |C,=K) 








L+t+T 
iE Ck, of tienen "i » (Pears =(£+t)') } 
StL+t+1 
= Cor 
A-¢ L+t+T 
+ Atk, | > (s+(t-k)*) - > (s+ (t-k)*) (1-@ (s-(L+t))) 
S=Lt+t+1 S=L+t+1 
J 
aac & (5) +(H-(E+t+L+J) ) 
s=1 
Beeee cL) (Kk, €) + B(K;t, €)n ,(k,t) - Clk;t, ¢t) (E.4) 
and 
Vai ee (tT, 6). (C.=K] 
= A(kjt, t)? y2(k, t) + B(kit, t)? vz(k, t) ay) 
+e Pe )Py ae) Veen) Ailkst,t)B(k;t,t) . 
Let 
C+1 
pty = omen Cc, CB (Ce,(t)-C,,(t, €) |Cp=k) . (E.6) 


k=1 
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then 
Var [Cop (t) - Cup (t, €)) 


C+1 
= Son" (k, t) Var [Cop t) — Gye (t, €) [Cp=K) 
k=1 


t+1 


Finally, the mean of the total cost advantage is 


Met, G) Tatinew, Ge Ge 


mt, €) 


+ 


s=1 


J 
coils [(-E=8-2) F a-ats))] 


J 
- Care Is a (s) + (H- ( t+t+L+J) |, 


S=1 


and variance of the total cost advantage is 
Var lO. (6) ae scr) 


= Var[Co.(t) -Cye(t, t)] 


+ Var[Ga,(t) -Eyy(t, €)] 


+ > n° (kK, t)(BICa.(£) — Cy. (ome) Gra Se Gate 
k=1 


(E.7) 


(E.8) 


(E.9) 


Since the other terms in the cost function are constant terms 


and we are assuming the failure time series and maintenance 


action time series are independent. 
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APPENDIX F 


DATA SAMPLE 


This Appendix contains a sample of the input data for both 
the decision aid developed in this thesis, the PASCAL program 
FCT.PAS and the decision aid currently in use, the ROI 
program. The data are for a radar transmitter for the F-14A. 
They are listed in the following tables. The ROI data is used 
directly in the program. The rest of the data sets are used to 
support the ROI data. 

Table 1 presents data used in the ROI program. All the 
values for the current subsystem except "$/Rpr" and "AVDLR 
Cost/Unit" are computed by taking average values of these 
measures for the last 24 months. They are computed by the ROI 
program automatically. The values for the improved subsystem, 
the cost for fix, and the lead time sections of the table are 
input by the analyst for the upgraded subsystem. The "other 
information" section of the table is computed by the ROI 
program. Table 2 presents flight hours and numbers of aircraft 
for each month ina 5 year period. Table 3 presents the number 
of BCM’s, the number of maintenance actions, the mean flight 
hours between failures, the mean flight hours between 
Maintenance actions, and the number of failures for each 


month. Table 4 contains the man-hours for failure, for 


ai 


organizational level maintenance action, and for intermediate 


level maintenance action for each month. 


A. ROI DATA 


TABLE 1. ROI DATA 





TEC/TMS: AFWA/F-14A 
WUC : 74A1500 - T1224/AWG9 RADAR TRANSMITTER 0 


Pie a a Current 
42.47 
20.60 
6.34 
9.01 
oe 
1082.64 
7.80 
1120.36 


1000000 
250000 
5000 
2500 
400000 
400000 
350000 


[Lead Time Mths 
18 
Start, dinsittilmofaekoses 3 


Total Lead Time 


Improved 
80.00 


20.60 
ee) 

« .00 
13 
1082.64 


1120035 


PR | th 
un fs 


| 
















Kits/Mth Installed 15 (Kits) 


Jeross over month | 8 
[Number of systems | 400 | 
ee ee = 1526 






MTBF (or MFHBF) - Mean time between failures 
MTBMA (or MFHBMA) - Mean time between maintenance actions 
"O" MH/MA - Organizational level maintenance Man 


-hours per maintenance action 

"I" MH/MA - Intermediate level maintenance Man- 
hours per maintenance action 

MH/F - Maintenance Man-hours per failure 

$/Rpr [Matl] - Cost of material (bit & piece cost) 
per repair 

# BCMs [1-8] /Mth - Average monthly number of BCMs 
categories 1-8 


AVDLR Cost/Unit - Unit cost per AV-DLR action 


i) 


B. FLIGHT HOURS AND NUMBER OF AIRCRAFT DATA 


TABLE 2. FLIGHT HOURS AND NUMBER OF AIRCRAFT 


BLT-HRS/ TOTsAear 
AFWA F-14A 


TOT-ACES 






oo fee 4241.0 


Toole 535 Ue 
7206. 








W JW 
NO fre 
nd a) 


199 E710 652657 


1991/04 
1991/08 


WwW 
th 
W 


326 


W 
bo 
WwW 


1990/07 
1990/06 
1991/05 
1991/04 402 
fisgo/o5 | CSS 
459 
472 | 
467 


O 
my 
ct 
o 

W 

WwW 

~ 


nay 
oO 
Ol 














[agsnoa | 10250.4] rn 
i Ss 
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ee sce | cc 
fisesya0 | ttt. as | 
eee |S tcons. 9] ew ase, 
eas 070.5 | 852 | 
eos | 88 | 
e970. | 27 | 
eesrou cero faa 
Bee] Ci | tt 
Mees] a 
fageeyog | ter. f 
fageeyoe | tag. | a 
meio ssf ta | 
fasessog | tag | a 
Bee tor] tt 
443 

ae  xut 


463 
458 

58 
458 
452 
452 
447 
441 
443 
447 
442 
448 
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1987/05 9445.1 
1987/04 9669.5 





C. BCM, MA, MFHBF, MFHBMA, VF DATA 


TABLE 3. BCM, MA, MFHBF, MFHBMA, VF 


BCM/MA/MFHBF /MFHBMA/VF 
AFWA F-14A 74A1500 - T1224/AWG9 RADAR TRANSMITTER 0 


MEHEMA 
faggasi2 | 3 | 20a] te. | 20.8 | al 
fagea/aa_| i | as 342s 
pasaayio | 9 |e eeaez | 36.2]| 18.0] Se 
luseayos | es | ass | 38.8) 20.0 | ney 
Pe ee ee 
fageayos | ss | a82| a2] 225 
lasgiyos | a0 aoa] a3 
50.3 
39.8 
42.9 
jagsovog | a] goa | 39.0] 20.8] 204 | 
jisgovoe | | soo | s9.6| 20.9 | | aes 
jissoo7 | is | asa] 0.7] to | 208 
fisgovos | sf aor] are] 3a] aa 
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fasgoo5 | a | aaa | ag | 22.2] 198 
ees | oy asa) tz.s | ars] a9 | 
jassgvoo | | 503 43.1 | 8.6] 247 
38.4| 18.7| 272 
fs 07 ie oil ae 
0. oa 
eee as eal alle 4.7 loss |e 20 
fee/00 7 oe 

es 

417 


3 
- 


1989/03 40.5 255 
1989/02 ee ss 210 
1989/01 407 41.8 236 
agee/ia | a] 70 asa | aaa] 196 | 
eee ei saz os] aa | 27 | 
essere ar | ars] tg.) 27 | 
faseayoo | a] ate | 3.7 | aaa | 237 | 
fageeyo7 | a] gov] a3 | ate | 2a 
aseeyos | 6 | aes | aaa] aa.o | 249 


1988/02 
1988/01 
sga7/a2 | at 50] az6] te] a7 








D. MH-FAILS, MH-ML1(2)-S, MH-ML1(2)-U0 DATA 


TABLE 4. MAN-HOURS FOR FAILURE AND MAINTENANCE ACTION 

















FWA F-14A 74A1500 - T1224/AWG9 RADAR TRANSMITTER 0 






Date MH-ML2-s | MH-ML2U | 
jigsi/i2 | e06.6| 0.6 | 1232.9) 0/0) Sema 
figsa/aa | 1952.9], 0.0 | 2420.4) ool) aaa 
jassa/1o | 2162.6] 0.0 |) 20ea.7/)_J0.0)| Saaaemey 
lissiyoa | 2642.5] 0.0] 2320.7] 0.0] 1566.1 
fassiyos | 3139.5] oo] 2760.5] 0.0] 2007.6 
fussaso7_ | 2227.9| 6.0 || Jaadia|) 762) Wasa 
jissiyos | 2470.8] 0.0] 2310.4] 0.0 | 1635.8 
jissiyoa | 2994.3] 0.0] 2522.0] 0.0] 1892.8 
jig91/o3 | 3170.5 o.0| 2690.0] 0.0] 1968.9. 
jassayo2z | 3953.6| 0.0] 3028.9]) 0.0 | 2eammay 
Jigoiyoa | aeea.a] 0.0] 3586.9] 0.0 | 3249.3 
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1990/11 | 2e71.1] 0.0 | 2576.0] oo] 1869.1 
fasso/io | 2616.4] 0.0] 2a5e.1| 0.0} 1695.5 
jisgovog | 2753.4] oo] 24sa.a] 0.0] 1724.3 
jisg0/os | 3603.1] 0.0] 2920.7] 0.0 | 2570.8 
1990/07 | 2734.1| 0.0] aera] 0.0] 1966.6 
jisgo/os | 2571.3] 6.0} 2107.8] 0.0 | 838.2 | 
jisgo/os | 2505.1] 2.4] 246.9 0.0] 1897.0 
fisso/oa | 2389.0] 0.0] 2344.9] 0.0] 1706.6 | 
jassovo3 | 3433.5] o.0| 3017.9] 0.0 | 2066.2 | 
pissoyo2 | 2653.0) oo} 2795.2] 0.0 | 1812.9 | 
jisgoso1_| 3405.1] 0.0] 3804.8] oo] 2154.9 
Bpeeyie | isc] Woo] ses] 0.0] a125.5) 
pises/ar 272.6] _o.0| 2577.6] 0-0 | 1282.5 | 
jisa9/o9 | 2705.6] 3.0] 2965.1] 0.0] 1704.9 | 
fasesvoe | 3302.6] 4.8] 3318.3] 4.8 | 2006.0 
piseeyo7 |] 2as7 a] ool) 2505.7} oo} 1497.7 
jises/os | 2906.1] 0.4] 2572.7] 0.0 | 1688.8 | 
jase9/os | 2881.9] 0.0] 2447.1] 0.0 | 1727.8 | 
jisesvoa | 3394.6] 0.0] 2895.9] 0.0 | 2255.0 
feeevos | 29523 oo) 2739.5] 0.0] 2015.4 
jases/o2 | 2550.0] 5.2] 2526.6] 0.0 | 143.8 
jisesvo1_| 3069.8] 12.6] 2582.5] 0.0] 1697.3 
j1gea/i2 | 2208.3 aa.a| 2387.2] 0.0] 1096.5 
jisea/ai | 2961.8] 0.0] 2695.2] 0.0] 1723.3 | 
jigee/1o | 3193.2] 0.0] 3634.9] 0.0] 2212.7 
jageesoo | 3009.2| 0.0] 2620.1] 0.0 | 1578.0 
jiseesos | 2817.8] 12.8] 2404.0] 0.0 | 1727.2 | 
jisee/o7 | 2632.7| oo] 2aez.2] 0.0 | 1645.9. 
jisse/o6 | 2192.8] 0.0] 2216.8] 0.0] 1503.5 
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198/05 | 3630.0] 0.0 | 3099.4 
19ge/o4 | 2689.5] 7.2| 2394.8 
198/03 | 2332.2] 7.0] 2626.2. 
1gse/oz | 3449.0] 3.5] 3207.2 
1988/01 16.6 
1987/12 
1987/11 


: : 2446.7 
1987/10 1942.0 


| 0.0 . 

| 0.0} 1469.5. 
| 0.0 | wenanal 
| 8.0 | asia 





| oo | 1464.0. 
| io. | 9200 
| 0.0 | 1e7eal 
| 0.0 | Team 
| oo | 1988.4 | 
| oo | 1780.3 | 
oo . 
rs 7 
pet . 
. 
pera 












Jane 

7.0 

Seo 
ao 
ao 
pee axo! 
19e7/o9 | 2433.2] 0.0] 2106.9 
19e7/08 | 3194.1] 0.0] 2449.1 
1987/07 | 3559.5] 0.0 2964.9 | 
1987/06 | 2647.4] oo] 2210.4 
1987/05) 295305) iba 
19e7/04 | 2054.5] 0.0. 
1987/03 | 2928.5| 0.0 
1987/02 
uge7/o1 | 2695.5] 0.0] 1832.6. 


The "O" level MH/MA are shown by its components MH-ML1-S 


2166.4 
2551.3 
221925 


1892.4 
173282 
192236 





and MH-ML1-U. These must be added together on a month by month 
basis to get the total "0" level MH/MA. Similarly the "I" 


level MH/MA is made up of MH-ML2-S and MH-ML2-U. 
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APPENDIX G 


PROGRAM OUTPUT 


A. OUTPUT FOR SIMULATED DATA 


1. Simulated data 


NO. of Feilures Meen ramber of Feilures 
(Maintenance Act ions ) of Maintenance Act ions 


Y-Axis scale 100 eee Press <Enter> to exit ss 


Figure 2. Simulated data. 





"Mean number of" "Mean number of" 

Failures Maintenance Actions 
219 403 442.669 

238 2545 45:0m67 9 

B21 .2038 4072159 
206.230 458.454 
230.087 426.491 
235.508 442.404 

213. Seg asa 
221.863 450.907 
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217 -8ou 454.586 


212,858 443.307 
214.141 4555779 
215. dom 457 . #50 
224.760 430,095 
229.2028 437.908 
237.5168 447.812 
220057. 451.614 
234.781 456.032 
203 2516 Assoc! 
Zee 418.431 
23652 0 441.045 
210.473 424.941 
Zee eeZ 449.765 
23. 0249)5 458.645 
202.616 0 443.813 
230.243 423.049 
223.042 437.962 
231.463 464.286 
236.364 466.401 
225 2209 £60 .691 
250 797 487.308 
261 793 456.121 
243 4863 480.219 
237016 A7a3792 
2567-560 494.986 
266.044 478.079 
284.469 488.819 
278.240 494.071 
280.863 Ao) oe 
263.594 497.664 
ZO 2.9 26 513.294 


2. Bayesian procedure 


Data (Read file: 0, Simulate: 1) 
Simu. para. (MUF, EtaF, Sig2F, CF ) 
Simu. para. (MuMA, EtaMA,Sig2MA, CMA) 
Reading filename 

NO. of systems & Use/Mon (NS.UPM) 
New MTBF & MTBMA (MTBFNew, MTBMANew) 
O MH/MA old & new (OMHOld,OMHNew) 
I MH/MA old & new (IMHOLd, IMHNew) 
MH/F old & new (FMHO1d, FMHNew) 
MATL Cost O & N (MATLO1d, MATLNew) 
NO. of BCMs O & N (BCMO1d, BCMNew) 
AVDLR Cost O & N (AVDLRO1d, AVDLRNew) 
Fix cost (CF) 

Lead Time & Kits/Mo (LeadTime, INSTL) 
Time start & horizon (ST,Hor) 
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1 
225 4 
440 a 


Oo ean ag Nae oa! 


400.00 
80.00 
6.34 
90a 
13 
1082.64 
7 FE6 
1120 336 
5400000 
60 

6 


144 25 
22a 22 


25n00 
20.60 
13. 
1082.64 
1120-35 


iS 
150 


Method (Likelihood: 
Replication (Rep) 


0, Bayesian: 1) : 1 


Cost advantaga 
of upgrade 


Assassins ee inty 
~ tine 6 40 





Y-Axis scale E+8 


Figure 3. 


#88 Press <Enter> to exit ss# 


Bayesian result for simulated data. 


Bayesian procedure 

The best time for subsystem upgrade 
and assessing uncertainty 

(from time 6 to time 40) 


"Time" "Best " " Cost " "Std.Dev." "Mean-2SD" "Mean+2SD" 
Index upgrade advantage Of “cost bound bound 
time of upgrade 

6 6 4.546E+07 8.266E+07 -1.199E+08 2.108E+08 

q 150 -2.614E+07 7.234E+07 -1.708E+08 1.1855+08 

8 8 1.305E+07 5.148E+07 -8.992E+07 1.160E+08 

9 150 -4.607E+06 4.636E+07 -9.734E+07 8.812E+07 

10 150 -1.937E+07 4.778E+07 -1.149E+08 7.619E+07 
apt 150 ~1.434E+07 3.520E+07 -8.475E+07 5 .607E+07 
a2 150 -9.970E+06 2.588E+07 -6.173E+07 4.179E+07 
des 150 -4.122E+06 2.527E+07 -5.467E+07 4.642E+07 
14 aA 7.006E+06 3.476E+07 -6.250E+07 7.652E+07 
15 15 3.878E+07 6.105E+07 -8.333E+07 1.609E+08 
16 16 1.056E+07 2.539E+07 -4.023E+07 6.134E+07 
17 17 1.968E+07 3.486E+07 -5.005E+07 8.940E+07 
18 150 -6.433E+07 9.237E+07 -2.491E+08 1.204E+08 
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19 150 -2.441E+07 3.863E+07 
20 20 1.567E+07 4.822E+07 
Dik TS0 ~1.315E+07 3.870E+07 
ae 150 -4.170E+06 2.271E+07 
23 26 6.650E+06 2.699E+07 
24 150 ~2.807E+07 6.101E+07 
25 150 ~1.630E+05 2.570E+07 
26 150 -6.507E+05 1.788E+07 
2a oe, 9.095E+06 2.684E+07 
28 28 1.845E+07 3.560E+07 
29 29 5.109E+06 1.667E+07 
30 30 6.800E+07 7.760E+07 
3. Sal 8 .249E+07 5 .462E+07 
32 a 3.344E+07 1 oO 
a 32 2.004E+07 1.003E+07 
34 34 2427 2E+00 8.962E+06 
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Figure 4. Likelihood result for simulated data. 
"Time" "Best " " Cost " "Std.Dev." "Mean-2SD" "Mean+2SD" 
Index upgrade advantage of cost bound bound 
time of upgrade 
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Figure 6. Bayesian result for real 


Data (Read file: 0, Simulate: 1) 
Simu. para. (MuF, EtaF, Sig2F, CF ) 
Simu. para. (MuMA, EtaMA, Sig2MA, CMA) 
Reading filename 

NO. of systems & Use/Mon (NS.UPM) 
New MTBF & MTBMA (MTBFNew, MTBMANew) 
O MH/MA old & new (OMHOld, OMHNew) 
I MH/MA old & new (IMHOLd, IMHNew) 
MH/F old & new (FMHO1l1d, FMHNew) 
MATL Cost O & N (MATLO1d, MATLNew) 
NO. of BCMs O & N (BCMOld, BCMNew) 


68 


40 30 


Mean = 
Maant+¢(-32SD 


60 


8 Press <Enter> to exit s## 


data. 


0 
0 


-b:\f15 corn 


400. 
80. 
Ge 

2 
132 
LOS2e 
Wir 


0 
0 
00 


00 
34 


Ot 


¥3 
64 
80 


25m 
205 
Ua 
OA 
oe 
.64 
.00 


13 
1082 





oOo 


00 
60 


AVDLR Cost O & N (AVDLRO1d, AVDLRNew) 
Fix cost (CF) 

Lead Time & Kits/Mo (LeadTime, INSTL) 
Time start & horizon (ST,Hor) 

Method (Likelihood: 0, Bayesian: 1) 
Replication (Rep) 


Bayesian procedure 

The best time for subsystem upgrade 
and assessing uncertainty 

(from time 6 to time 60) 


fsieme" “Best “" "“ Cost ** ™Std.Dev." 
Index upgrade advantage O£ cost 
time of upgrade 

6 180 -3.031E+08 3.154E+08 
7 7 1.583E+08 3.020E+08 
8 8 4.605E+07 1.692E+08 
9 180 -~1.832E+08 2.988E+08 
10 180 -2.162E+08 2.571E+08 
ale 180 -8.537E+07 1.079E+08 
Ale 2 6.514E+08 4.816E+08 
13 13 4.591E+08 3.398E+08 
14 14 7.408E+07 1.270E+08 
15 180 -8.311E+07 2.696E+08 
16 180 ~4.158E+07 1.589E+08 
11% 7 3.493E+07 1.510E+08 
18 180 -1.397E+08 3.003E+08 
19 19 5.106E+07 1.915E+08 
20 180 -1.503E+07 1.259E+08 
21 2. 1.374E+07 1.218E+08 
2 Ze 2.720E+08 3.470E+08 
23 23 9 .428E+07 1.649E+08 
24 Za 3.181E+07 8.776E+07 
25 25 3.559E+07 9.217E+07 
26 216 1.502E+07 726528707 
a7 24 3.652E+07 9.307E+07 
28 28 7.122E+07 1.333E+08 
29 29 2.423E+07 4.937E+07 
30 30 3.123E+07 6.212E+07 
ool shih 2.130E+07 3.823E+07 
eZ a2 3 .330E+07 6.929E+07 
eo 33 2.517E+07 2.952E+07 
34 34 2.187E+07 2.656E+07 
35 35 2.335E+07 2.654E+07 
36 36 2.294E+07 2.492E+07 
a7 37 3.770E+07 5.581E+07 
38 38 2.821E+07 1.798E+07 
39 39 2.407E+07 2.013E+07 
40 40 1.301E+07 4.390E+07 
41 41 7.387E+06 4.962E+07 


69 


1L1i20 356 
5400000 
60 

6 

al 

100 


“Mean-2SD" 


=o 
-4 
a 
=e 
=i) 


=3 
= 
=H 
=0 
=o 


= 7 
= 
a 
= 


el 


aoa 
2 
—« 
= 
-@ 
=5 


=3 
=5 
ae 


bound 


.338E+08 
.456E+08 
.923E+08 
.808E+08 
.304E+08 
-3. 


013E+08 


.118E+08 
.206E+08 
.798E+08 
.223E+08 
.595E+08 
-2. 
.404E+08 
.320E+08 
.668E+08 
.299E+08 
-4, 
-2. 


670E+08 


220E+08 
355E+08 


.437E+08 
-1. 
. 380E+08 
.496E+08 
.953E+08 
-450E+07 
.301E+07 
.516E+07 
-i. 


487E+08 


053E+08 


.388E+07 
.126E+07 
.973E+07 
~2. 
.392E+07 
.757E+06 
.620E+07 
-479E+07 
. 186E+07 


691E+07 


B20. 36 


WS) 
180 


"Mean+2SD" 


PRPAHDRPITITWTORPWOPRPPWNHRNHYNBPUWONNFKBWNHKPWRPPPRPN BW IW 


bound 


.276E+08 
.622E+08 
.844E+08 
.145E+08 
-980E+08 
.305E+08 
.615E+09 
.139E+09 
.280E+08 
.561E+08 
.763E+08 
.369E+08 
.610E+08 
.341E+08 
.368E+08 
.574E+08 
.660E+08 
.240E+08 
.073E+08 
.199E+08 
.681E+08 
.227E+08 
.378E+08 
.230E+08 
.555E+08 
.776E+07 
.719E+08 
.421E+07 
-499E+07 
.643E+07 
.279E+07 
.493E+08 
-417E+07 
.433E+07 
.008E+08 
.066E+08 


42 42 6.332E+06 4.417E+07 -8.201E+07 9 .467E+07 
43 43 1.468E+07 3.849E+07 -6.231E+07 9.167E+07 
44 44 1.615E+07 3.932E+07 -6.249E+07 9.479E+07 
45 aS 6.863E+06 3.979E+07 -7.272E+07 8.645E+07 
46 180 ~3.456E+06 6.481E+07 -1.331E+08 1.262E+08 
rl) A 1.015E+07 3.697E+07 -6.379E+07 8.409E+07 
48 48 9.418E+06 3.168E+07 -5.394E+07 7.278E+07 
49 49 1.316E+07 4.153E+07 -6.990E+07 9 .623E+07 
50 180 -1.373E+06 5.540E+07 -1.122E+08 1.094E+08 
5 Si 5 .869E+06 2.976E+07 -5.365E+07 6 .539E+07 
52 5 8 .059E+06 3.121E+07 -5.437E+07 7.049E+07 
55 55 1.619E+06 3.606E+07 -7.051E+07 7.375E+07 
54 54 6.219E+07 1.549E+08 -2.475E+08 3.719E+08 
SS a5 1.092E+08 1.535E+08 -1.979E+08 4.163E+08 
56 56 5.206E+07 7.508E+07 -9.810E+07 2.022E+08 
57 57, 2.849E+07 4.166E+07 -5.483E+07 1.118E+08 
58 58 2.637E+07 3.311E+07 -3.986E+07 9.260E+07 
59 59 3 .835E+07 3.285E+07 -2.735E+07 1.040E+08 
60 60 9.829E+06 7.842E+06 -5.855E+06 2.551E+07 


3. Maximum likelihood procedure 


Cost advantaga Asaaasing uncertainty 
cf, 6 to 60 °* 


from tine 


Maan —__ 
Heant+¢(->32S0 ——— 





V-Axis scale E+9 


Figure 7. 


+t Press <Enter> to exit sss 


Likelihood result for real data. 


7/0 


Data (Read file: 0, Simulate: 1) 
Simu. para. (MuF, EtaF, Sig2F, CF ) 
Simu. para. (MuMA, EtaMA, Sig2MA, CMA) 
Reading filename 


NO. of systems & Use/Mon (NS.UPM) 

New MTBF & MTBMA (MTBFNew, MTBMANew) 
O MH/MA old & new (OMHO1ld, OMHNew) 
I MH/MA old & new (IMHOLd, IMHNew) 
MH/F old & new (FMHO1ld, FMHNew) 
MATL Cost O & N (MATLO1d, MATLNew) 
NO. of BCMs O & N (BCMOl1d, BCMNew) 


AVDLR Cost O & N (AVDLRO1d, AVDLRNew) 
Fix cost (CF) 

Lead Time & Kits/Mo (LeadTime, INSTL) 
Time start & horizon (ST,Hor) 
Method (Likelihood: 0, Bayesian: 
Replication (Rep) 


1) 


Maximum likelihood procedure 

The best time for subsystem upgrade 
and assessing uncertainty 

(from time 6 to time 60) 


0 

0 0 

; 0 0 
on \ £5 oven 
400.00 

B02 00 

6.34 

9 30H. 

i Sreies 

1082 ed 

7 280 
i206 
5400000 

60 

6 

0 

100 


“"Mean-2SD" 
bound 


eiime® “Bestetee"™ Cost ™ "Std.Dev." 
Index upgrade advantage Of Cosi: 
time of upgrade 

6 180 -5.224E+08 3.242E+08 
7 7 4.100E+08 3.364E+08 
8 8 2.383E+08 2.984E+08 
9 180 -4.743E+08 3.667E+08 
10 180 -4.726E+08 2.750E+08 
igh 180 -3 .022E+08 2.540E+08 
2 ee 9.329E+08 4.386E+08 
le ie 7.144E+08 4.120E+08 
14 14 3.542E+08 3.189E+08 
85 180 -4.399E+08 4.887E+08 
16 180 -2.532E+08 2.931E+08 
ley a, 1.832E+08 3.368E+08 
18 180 -4.825E+08 4.317E+08 
19 19 3.658E+08 4.353E+08 
20 180 -8.019E+07 3.485E+08 
Dei 21 1.228E+08 3.198E+08 
22 22 5.815E+08 3.718E+08 
23 23 3.945E+08 2.859E+08 
24 24 1.514E+08 2.211E+08 
25 25 1.825E+08 2.479E+08 
26 26 1.375E+08 2.518E+08 
27 27 1.070E+08 1.742E+08 
28 28 2.338E+08 2.335E+08 
29 29 1.203E+08 2.037E+08 
30 30 9.947E+07 1.792E+08 


We 


al 
=e 
3 
=i 
aoe 


=2 
=a 


-8., 
.905E+08 
-1. 
-5. 
.772E+08 
.167E+08 


-4 


= 
aS 


=e 
.773E+08 
. 908E+08 
.133E+08 
.662E+08 
.414E+08 
.333E+08 
.870E+08 
.588E+08 


=. 
2 
ao 
=3 
=e 
=2 
ae 
ma, 


.171E+09 
.628E+08 
.585E+08 
.208E+09 
.023E+09 
-8. 

Bs 
-1. 
.836E+08 
.417E+09 


102E+08 
578E+07 
097E+08 
394E+08 
346E+09 
049E+08 


621E+08 


Ze. 
ZU. 


00 
60 


13 
1082 


pales 
64 
120235 


15 
180 


"Mean+2SD" 


Em UTA DMHDUWOPAATIHDRPWOWUWORFPERPNATIN OFF 


bound 


-.260E+08 
.083E+09 
.352E+08 
.591E+08 
.748E+07 
.0O58E+08 
.810E+09 
~538E+09 
.919E+08 
.376E+08 
.3I30E+08 
-569E+08 
.808E+08 
.236E+09 
.168E+08 
.624E+08 
-325E+09 
.664E+08 
.936E+08 
. 782E+08 
.411E+08 
.555E+08 
.009E+08 
.277E+08 
.578E+08 


WIDRPrPWWWATNHONRPWWPNHKRKRWUBPNPHDOUP UN UI 


.994E+07 
.364E+07 
-.474E+07 
.988E+07 
.438E+07 
- 332E+07 
.361E+07 
-272E+07 
. 866E+07 
.781E+07 
.147E+07 
.562E+07 
-570E+07 
-500E+07 
~252E+07 
.480E+08 
.083E+07 
1126407 
-914E+07 
.607E+07 
.750E+07 
.305E+06 
.637E+07 
.591E+08 
.088E+08 
.464E+08 
.027E+08 
.427E+07 
.846E+07 
.585E+07 


MOYYORPNPRPRPRPRPRPRPRPPOPRPPRPEFPYIOWOOMPRPIYIOP 


.046E+08 
.879E+07 
-949E+07 
.618E+08 
.703E+07 
.323E+07 
.020E+07 
.133E+07 
.017E+08 
.511E+08 
.126E+08 
.164E+08 
.229E+07 
.263E+08 
.419E+08 
.886E+08 
.419E+08 
.360E+08 
.288E+08 
.104E+08 
.423E+08 
.204E+08 
.318E+08 
.313E+08 
.963E+08 
.983E+07 
.619E+07 
.684E+07 
.815E+07 
.340E+07 


Te 


.494E+08 
.139E+08 
.042E+08 
. 638E+08 
-297E+08 
.331E+08 
.168E+08 
.994E+07 
. 748E+08 
.544E+08 
.736E+08 
.972E+08 
.389E+08 
.076E+08 
.612E+08 
.251E+08 
.531E+08 
.408E+08 
.385E+08 
.947E+08 
.571E+08 
.336E+08 
.272E+08 
*,034E+08 
. 384E+07 
.321E+07 
-965E+07 
-942E+07 
.784E+07 


094E+07 


PMN NWATOMOWNWNHNWWNHOWNHDNNNWNHRPNYNNNWNN ND 


.692E+08 
.412E+08 
.137E+08 
.835E+08 
.184E+08 
-398E+08 
-440E+08 
.854E+08 
.321E+08 
-500E+08 
. 766E+08 
-684E+08 
.303E+08 
-976E+08 
.062E+08 
.291E+08 
.147E+08 
.031E+08 
. 768E+08 
.468E+08 
.121E+08 
.482E+08 
.0O00E+08 
.217E+08 
.015E+08 
.461E+08 
.551E+08 
.280E+08 
-548E+08 
-.426E+08 


APPENDIX H 


UPGRADE.PAS PROGRAM 


A. USER DOCUMENTATION 
This documentation contains the information concerning the 
utilization, input data, and results of the UPGRADE. PAS 
program. 
1. Utilization 
This decision aid program (UPGRADE.PAS) is developed 
to enhance the Economic Analysis program (ROI) in the 
Automated Management Indicator System (AMIS). It uses 
statistical procedures to estimate the time of onset of 
subsystem degradation and the magnitude and evolution of the 
degradation over time. These estimates are used to compute the 
estimated cost of remaining with the current subsystem. A 
comparison of this cost with the cost of investing in the 
upgraded subsystem can suggest to the user a most economical 
time to upgrade this particular subsystem. The analysis 
includes assessment of the uncertainty in the estimated costs. 
This program is based on the model, which postulates 
that there may be a linear trend in the mean number of 
failures (maintenance actions) per time period for the current 


subsystem. It uses a statistical method to detect and quantify 
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the trend and combines these with the other values to estimate 
the future costs for the subsystem. 

The program is written in TURBO PASCAL for use on a 
personal computer, not for the mainframe. It also includes 
screen graphics. The default screen graphics is VGA. The user 
Should adjust the value of "GraphDriver", "GraphMode", and the 
path of "InitGraph" in procedures "GenDataGraph" and 
"GenGraph" before using the program. The graphical output 
cannot be printed by pressing <Print Screen> on keyboard. The 
user needs other software (Colorix etc.) to get a hard copy of 
the graphs. 

2. Menu and Input Data 

When the program is run, the menu will appear on the 
screen and wait for further adjustments to default values of 
the parameters. Permanent changes to parameter values should 
be done in the procedure "DefaultParameter" of the source 
code. All input data was chosen to be similar to those used in 
the ROI program, so the values of the parameters can be 
obtained from other existing data bases (NALDA etc.). The main 
difference between the UPGRADE.PAS and ROI programs is the 
estimation of the mean number of failures (respectively 
maintenance actions) for the current subsystem. In UPGRADE. PAS 
we use the whole collected data array to detect and quantify 


the trend. The resulting estimates are used to estimate the 
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future costs. The ROI program simply uses the average values 


of the measures of performance for the past 24 months data. 


A listing of the menu, excluding default values for 


the parameters, is 


0 


1 


10 
Bal. 
a2 
3 
14 
15 
16 
i) 


18 


Run the model. 

Data (Read file: 0, Simulate: 1) 
Simu. para. (MeuF, EtaF, Sig2F, CF ) 
Simu. para. (MeuMA, EtaMA, Sig2MA,CMA) 
Reading filename 

NO. of systems & Use/Mon (NS.UPM) 
New MTBF & MTBMA (MTBFNew, MTBMANew) 
O MH/MA old & new (OMHO1d, OMHNew) 
I MH/MA old & new (IMHOLd, IMHNew) 
MH/F old & new (FMHO1d, FMHNew) 
MATL Cost O & N (MATLO1d, MATLNew) 
NO. of BCMs O & N (BCMO1d, BCMNew) 
AVDLR Cost O & N (AVDLRO1d, AVDLRNew) 
hax Cost» (CE) 

Lead time & Kits/Mo (LeadTime, INSTL) 
Timemsibantuc&ehorizone(ST,Hor) 

Method (Likelihood: 0, Bayesian: 1) 
Replication (Rep) 

Reset to default values. 


Items 1 through 4 ask the user to choose the data 


array of mean number of failures (respectively maintenance 


actions) for the current subsystem (item 1). If the user 


de 


elects to simulate these two arrays, then he/she should input 
the values of the four model parameters for each array; these 
are the mean value before the trend, the slope of the trend, 
the standard deviation of the data array, and the occurrence 
time of the trend (items 2 and 3). Item 2 requests the 
parameters used for the number of failures for the simulation; 
item 3 requests the parameters used for the number of 
maintenance actions for the simulation. If the user decides to 
analyze outside data, a ASCII file is needed containing the 
MFHBF and MFHBMA arrays in the NALDA data base (item 4). The 
format of the file should have two columns. The first column 
is MFHBF; the second column is MFHBMA (see columns 4 and 5 in 
Table 3). The UPGRADE.PAS program will transform these two 
arrays to the mean number of failures and maintenance actions. 

Items 5 through 12 can be obtained from the ROI 
program in the "Current System" and "Improved System" 
sections. The user inputs the same values for the identified 
items. The UPGRADE.PAS program is designed to read the 
pairwise parameters (both systems) together. The parameters 
have the same heading followed by "Old" for the current system 
parameters; those followed by "New" are the improved system 
parameters. 

Items 13 and 14 can also be obtained from the ROI 
program in the "Cost of Fix" and "Schedule for Fix" section. 
The items "Cost/Kit" and "Cost for INSTL/Kit" should be 


multiplied by number of systems (item 5) and added to the 
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other items in the cost of fix section to get "CF" (item 13). 
To compute "LeadTime", the user needs to add all the items 
together except the last item "Kits/MO INSTL" in the "Schedule 
for Fix" section. The value of "INSTL" is identical to that of 
"Kits/MO INSTL" in the ROI program (both at item 14). 

Items 15 to 17 ask the user to choose the desired time 
in the time series to start the estimation of the best 
upgrading policy and mean cost advantage, and the expected 
time horizon for the subsystem (item 15). Then it provides two 
methods (the Maximum Likelihood and Bayesian) to estimate 
costs. The Bayesian method requires less computation than the 
Maximum Likelihood method (item 16). To only obtain the best 
policies for each time for the chosen procedure, set the 
parameter Rep=0. If Rep is a positive integer, then 
assessments of uncertainty for the best policy will also be 
given. If the Bayesian procedure is chosen, the assessment of 
uncertainty uses moments of the posterior distribution. If the 
maximum likelihood is used, Rep is equal to the number of 
replications for the bootstrap estimates of uncertainty (item 
17). If the Bayesian procedure is chosen the time to start the 
calculation, ST=6. 

The values of "MA/ML1(2) Ratio" and "ML1(2) Rate" can 
only be changed in the source code. They are stored in the 


constant declaration of the main program. 


wey 


3. Results 

When the program executes (choose item 0), three 
sections of results will appear. First, the graphs of the mean 
number of failures and mean number of maintenance actions will 
appear on the screen. Then, the most economical time for 
upgrading and estimated cost advantage of upgrade and two 
standard deviation bounds will appear after each time index 
(from the time chosen to start the calculation, ST, to the end 
of the time series, ET). Finally, the graphs of the mean cost 
advantage and two standard deviation bounds will appear on the 
screen. After looking the graph, user may press <Enter> to 
return the original screen. 

The results, excluding the graphics, will be located 
in the F.OUT file. Simu.Data file contains the simulated data 
(mean time between failure and mean time between maintenance 


action) for the purpose of reuse. 


B. SOURCE CODE 


{SM 36384,0, 655360} 
PROGRAM UPGRADE; 
Uses Dos, Crt, Graph; 


{ This program is developed to estimate the time of onset of 
a given (current) subsystem degradation and the magnitude of 
the degradation. These estimates are then used to estimate the 
cost of remaining with the current subsystem for the remaining 
time horizon, the life of the parent system. We compare this 
cost with the cost of investing in the upgraded (improved) 
subsystem to obtain a best time to invest in the upgraded 
subsystem. 

Three options are given in this program: 
1. Data acquisition is read from a data file (set Data = 0) 
or simulated from the program itself by choosing parameters 
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(Bacaw=. 1). 

2. Estimation uses the Maximum Likelihood procedure (set Meth 
= 0) or the Bayesian procedure (Meth = 1). 

3. Assess uncertainty or not. Set (Rep = 0) when the answer 
1s no. Otherwise the value of Rep represents the number of 
bootstrap replication. The Bayesian procedure will 
automatically assess uncertainty when Rep > 0 and Meth = 1. 


All input data should be copied from the ROI program 
except the MFHBF and MFHBMA arrays. These arrays should be 
contained in input file. They are obtained from the NALDA data 
base. 

The values of "MA/ML1(2) Ratio" and "ML1(2) Rate" can only 
be changed in the source code. They are stored in the constant 
declaration of the main program. The other parameters can be 
adjusted in the menu given on the computer screen. 

The length of "IntVec" and "RealVec" should be larger than 
or equal to the value of "ET+1". 

The length of "RealVecl" should be larger than or equal to 
the value of "Rep". 


When this program running, the graph of the mean number of 
failures and mean number of maintenance actions will appear on 
the screen. It is followed by the estimated best time for 
upgrading the subsystem that will be printed after each 
decision time index. 

If the user chooses to assess the uncertainty of 
estimation, then the cost advantage of the best upgrade policy 
and two standard deviation bounds will be printed on the 
screen following the best upgrading time. The user can view 
the graph and press <Enter> to leave graphical screen. 


The results, excluding the graphics, are located in F.OUT 
(Output) file. This file also contains the desired parameters 
(menu) for the computation. 

Simu.Data (Output2) contains the simulated data for the 


purpose of reuse if Data = 1. } 
const MAML1Ratio = 0.92; 

MAML2Ratio = 0.47; 

ML1Rate = 15.226 

ML2Rate = 18.355 


type IntVec array[1..121] of integer; 


RealVec = array[1..121] of real; 
RealVecl = array[1..1000] of real; 
StrvVec = array[0..18] of string; 
var CTVec > Inevec; 
NOFVec, NOMAVec : RealVec; 


ne 


SimuFVec, SimuMAVec RealVec; 
MuF, EtaF, SigmaSqrF RealVec; 
MuMA, EtaMA, SigmaSqrMA RealVec; 
VSqrF, KSquFk  skmecanh RealVec; 
VSqrMA, RSqrMA, RhoSqrMA RealVec; 
PiStarF, PiStarMA RealVec; 
PiDistF, PiDistMA RealVec; 
Mean, LB, UB RealVec; 
Gain RealVec1; 
Mission Strvec; 
RanSeed Longant, 
Data, Meth integer; 
ST, ED Hor integer; 
LeadTime, INSTL integer; 
Rep, R, Time, CT integer; 
SetCF, SetCF1l, HatCF integer; 
SetCMA, SetCMA1, HatCMA integer; 
MTBFNew, MTBMANew real; 
NOFNew, NOMANew real; 
OMHOl1ld, OMHNew real; 
IMHOld, IMHNew real; 
FMHOld, FMHNew real; 
MATLO1d, MATLNew real; 
BCMOld, BCMNew real; 
AVDLROld, AVDLRNew real; 
COMA, COF, COAD real; 
CNMA, CNF, CNAD real; 
CF, NS, UPM real; 
SetMuF, SetEtaF real; 
Set SigmaSqrF real; 
SetMuFl, SetEtaFl real; 
SetSigmaSqrF1l real; 
HatMuF, HatEtaF real; 
Hat SigmaSqrF real; 
SetMuMA, SetEtaMA real; 
Set SigmaSqrMA real; 
SetMuMA1, SetEtaMA1 real; 
Set SigmaSqrMA1 real; 
HatMuMA, HatEtaMA real; 
HatSigmaSqrMA real; 
Max, Min, Total real; 
Infile string? 
OuEDUIE, OuEpUEZ text; 
Input text ; 
{+++PART1 SET MENU & INPUT DATA4+4+44+44+4+4+4+44+4+4+4+4+4++4++4++4 } 
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{ These procedures set up the menu on the screen and input 
values for the parameters. It includes: 

Default Parameter, 

SetUnitCost, 

SetMission, 

SetParameter, 

WriteParameter. } 


Procedure Default Parameter; 


begin 
Data =) 0? 
Set MuF1 wo oe 
SetEtaFl = a0 
SetSigmaSqrFl = 225 
SetCFl1 come 1 ()'- 
SetMuMA1 = 440; 
SetEtaMAlL eer 
SetSigmaSqrMA1 := 400; 
SetCMA1 = 10% 
NS = 400; 
UPM =_ 2 Se 
Infile = 'a:\£15 yemen’; 

{MTBFOld, MTBMAO1d} 

OMHO1d = 6547 
IMHOld =e Uae. 
FMHOld — is ach sie 
MATLO1d = 1082.64; 
BCMOld =e 7 10r- 
AVDLRO1d = 1120.36; 
MTBFNew =o Ue O- 
MTBMANew = 208 60- 
OMHNew 9.0 ly 
IMHNew =D 20 las 
FMHNew ae oes eo 
MATLNew = 1082m64; 
BCMNew = Oe 
AVDLRNew =a? 0. 36> 
CF = 5400000; 
LeadTime = 60; 
INSTL eels 
ST =) iG 
Hor = ASOT 
Meth ele 
Rep = Oe 
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end; 


LHEFEHEEEEEEEHHEE HEHEHE HEHEHE EEE EEE EHH FHF +H $F +ttttget tet teet+t++44 i 
{ This procedure uses the input data to compute the six unit 
costs: } 


Procedure SetUnitCost; 
var OMHCostOld, IMHCostOld, FMHCostOld : real; 
OMHCostNew, IMHCostNew, FMHCostNew : real; 


begin 
NOFNew 
NOMANew 


NS*UPM/MTBFNew; 
NS *UPM/MTBMANew ; 


nou 


OMHCostOld OMHO1d*MAML1Ratio*ML1Rate; 
IMHCostOld IMHO1d*MAML2Ratio*ML2Rate; 
FMHCostOld FMHOl1d*ML2Rate; 

COMA OMHCostOld+IMHCostOld; 

COE FMHCostOl1ld+MATLO1d; 


COAD BCMO1d*AVDLRO1d; 


i "ou 


OMHCostNew OMHNew*MAML1Ratio*ML1Rate; 
IMHCostNew IMHNew*MAML2Ratio*ML2Rate; 
FMHCostNew := FMHNew*ML2Rate; 
CNMA := OMHCostNew+IMHCostNew; 
CNF FMHCostNew+MATLNew; 
CNAD BCMNew* AVDLRNew; 

end; 


LHAFEEHEEFHEEEFEFEEEEEEFEEEEEFHEEE HEE EE EF EFF + 4¢4¢444¢4+44+4+4+4+4+44+444 } 
Procedure SetMission; 


begin 
Mission[0] 
Mission([1] 
Mission [2] 
Mission [3] 
Mission[4] 
Mission[5] 
Mission[6] 
Mission[7] 
Mission[8] 
Mission[9] 
Mission[10] 
Mission[11] 
Mission([12] 
Mission[13] 
Mission([14] 
Mission[15] 
Mission[16] 


~ 


Run the model.’ 

Data (Read file: 0, Simulate: 1) 
Simu. para. (MuF, EtaF, Sig2F, CF ) 
Simu. para. (MuMA, EtaMA,Sig2MA, CMA) 
Reading filename 

NO. of systems & Use/Mon (NS.UPM) 
New MTBF & MTBMA (MTBFNew, MTBMANew) 
O MH/MA old & new (OMHO1ld, OMHNew) 
I MH/MA old & new (IMHOLd, IMHNew) 
MH/F old & new (FMHO1ld, FMHNew) 
MATL Cost O & N (MATLO1d, MATLNew) 
NO. of BCMs O & N (BCMOld, BCMNew) 
AVDLR Cost O &N (AVDLRO1ld, AVDLRNew) 
Fix cost (CF) 

Lead Time & Kits/Mo (LeadTime, INSTL) 
Time start & horizon (ST,Hor) 
Method (Likelihood: 0, Bayesian: 1) 


~ - 


a) 


Hou d uw dm bd db db kb td bem bb egy 
eo ee eo ee @ 
~ ~ - ~ ~ ~ ~ ~ ~ -~ ~ ~ -~ ~ ~ - 
™e oe Me Me Se Be NO Be Ne Be Me We We BO Be WE 


~- = ~ =~ -~ ~ ~ ~ ~ ~ ~ - = 


a2 


‘ Replication (Rep) oe; 
‘ Reset to default values.’; 


Mission([(17]} 
Mission[18] 
end; 


LEEFEEEEEEEEE EEE EEE EEE EEE EEE HEHEHE EHH E+E +4 4444444 4+4+4+4+4+4+4+4+4+4+4+4+4+4} 
{ This procedure allows the user to choose the desired 
values for each parameter. 

{ If Bayesian method (Meth = 1) is chosen, then starting 
time will set to 6 automatically. } 


Procedure SetParameter; 
var Choise, D, M : integer; 


{HEHEHE HEHEHE EHH HHH +4 et teeeetteeett++tses+ i} 


Procedure PrintMenu; 
var Choise : integer; 


begin 

writeln(’If you want to change the values of any 
parameter, '); 

writeln(’enter the number from nemu, ’); 

writeln(’or enter 0 for running the model.’); 

fo~ Choise <2°0Mto 18 do 

writeln(Choise:2,’ ‘',Mission[Choise]); 

end; 


{FAEEHHEEEE EHH HEEEE HEHE +4444 ++4++4+++4+4+4+4+4} 
Procedure ChooselItem; 


begin 
repeat 
readln(Choise) ; 
if (Choise < 0) or (Choise > 18) then 
GotoxM@y23)% 
until (Choise >= 0) and (Choise <= 18); 


if (Choise <> 0) then 
write(’Enter’,Mission[Choise]} ) 


else 

if (Data = 1) then 

begin 
writeln(’How many data do you want 

in the arrays (ET)? ‘'); 

write(’ (Beware ET <= Mor-LeadTime) ‘); 
readln (ET) ; 

end; 

end; 


LHEFEEEEEHHEE EEE FE FEE HE FE ++ +++ +++ +++++} 
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Procedure SetInitial; 


begin 

if Data = 0 then 

begin 
SetMuF 
SetEtaF 
SetSigmaSqrF 
SetCF 
Set MuMA 
SetEtaMA : 
SetSigmaSqrMA : 
SetCMA 

end 

else 

begin 
set MuF 
SetEtaF 
SetSigmaSqrF 
SetCFr 
SetMuMA 
SetEtaMA 
SetSigmaSqrMaA : 
SetCMA 

end; 

if Meth = 1 then 
Sl >-=o0, 


me ™e Me Me MO MO BO MS 


tow ot tw Ww mt tt 
oOo0oo0o0o0o0c00°0o 


SetMuF1; 
SetEtaFl; 
SetSigmaSqrF1; 
SetCF1; 
SetMuMAI]1 ; 
SetEtaMA1; 
SetSigmaSqrMAl ; 
SetCMAI] ; 


iow tt WoW om Ww tt 


end; 
{+E +EEEEEFHEE E+ HH t+ tetettttt+et+tt++t++)} 
procedure PrintParameter; 


begin 

GotoxXy (45,5); write(Data:10) ; 

GotoXy(45,6); write(SetMuF:5:0,SetEtaF:5:0, 
SetSigmaSqrF:5:0,SetCF:5); 

GotoxY(45,7); write(SetMuMA:5:0,SetEtaMA:5:0, 
SetSigmaSqrMA:5:0,SetCMA:5) ; 

GotoXYy (45,8); write(Infile:10) ; 

GotoxXY (45,9); write(NS:10:2,UPM:10:2) ; 

GotoxY (45,10); write (MTBFNew:10:2,MTBMANew:10:2); 

GotoxXY(45,11); write (OMHO1d:10:2,OMHNew:10:2) ; 

GotoXY (45,12); write(IMHO1d:10:2, IMHNew:10:2); 

GotoxY (45,13); write(FMHO1d:10:2,FMHNew:10:2); 

GotoxY (45,14); write (MATLO1d:10:2,MATLNew:10:2) ; 

GotoxXY (45,15); write(BCMO1d:10:2,BCMNew:10:2) ; 

GotoxXY (45,16); write (AVDLRO1d:10:2,AVDLRNew:10:2) ; 

GotoXY (45,17); write(Cr-10- on. 

GotoxYy (45,18); write(LeadTime:10,INSTL:10) ; 

GotoxY (45,19): write(ST:10,Hor:10) ; 
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GotoxXY (45,20); write(Meth:10); 
GotoXY(45,21); write(Rep:10); 
GoeeKY(1,23); write? ‘); 


end; 


{HHH +44 4444+ +++ +4+4¢4¢+4¢+4+4+4+4+4+4+4+4+4+4+4+4+4+4+44 } 


begin {Procedure SetParameter} 


repeat 
Cirsen: 
PrintMenu; 
SetInitial; 
PrintParameter; 
ChooselItem; 
case Choise of 
0 Goto X Yalen24) 
a begin 
readlin(D); 
if (Db = Oi or (D =a) then 
Daca o= 1): 
end; 
2 readln(SetMuF, SetEtaF, SetSigmaSqrF,SetCF) ; 
3 readln(SetMuMA, SetEtaMA, SetSigmaSqrMA, SetCMA) ; 
4 readln(Infile) ; 
5 readln(NS,UPM) ; 
6 readln (MTBFNew, MTBMANew) ; 
7 readin (OMHO1d, OMHNew) ; 
8 readln(IMHOld, IMHNew) ; 
9 readln (FMHOld, FMHNew) ; 
vO readln (MATLO1d, MATLNew) ; 
eal readln (BCMOld, BCMNew) ; 
2 readln (AVDLRO1d, AVDLRNew) ; 
aS} readln(CF); 
14 readln(LeadTime, INSTL) ; 
cS readln(ST,Hor); 
iG begin 
readin(M); 
if (M = 0) or (M = 1) then 
Meth := M; 
end; 
7 readln (Rep) ; 
18 DefaultParameter; 
end; 


until (Choise = 0); 


end; 


{Procdure SetParameter} 


LEFEEEFEFEEHEFE FE FE FEE FTE EE EE FEE EEE EEFE FETE HEHE +t +t t+t+t+t+t++++ +} 


Procedure WriteParameter; 


begin 


So 


writeln(Output,Mission[1]:10,Data:10); 

writeln(Output,Mission[2]:10,SetMuF:5:0,SetEtaF:5:0, 
SetSigmaSqrF:5:0,SetCF:5); 

writeln(Output,Mission[3]:10,SetMuMA:5:0,SetEtaMA:5:0, 
SetSigmaSqrMA:5:0,SetCMA:5) ; 


writeln(Output,Mission[4] 
writeln(Output,Mission[5] 
writeln(Output,Mission[6] 
writeln(Output,Mission[7] 
writeln(Output,Mission[8] 
writeln(Output,Mission[9] 


: 10OPi nla: 
:10,NS:10:270PM- Toe 2) 
:10,MTBFNew:10:2,MTBMANew: 
:10,OMHO1d:10:2,OMHNew:10 
2:10, IMHO1d:10:2, IMHNew:10: 
:10,FMHO1d:10:2, FMHNew:10: 


10 2a 
vane 
25 
2)% 


writeln(Output,Mission[{10] :10,MATLO1d:10:2,MATLNew:10:2) ; 
writeln(Output,Mission[{11] :10,BCMO1d:10:2, BCMNew:10:2); 
writeln(Output,Mission[12] :10, AVDLRO1d:10:2,AVDLRNew:10:2) ; 


writeln(Output,Mission[13] 
writeln(Output,Mission[14] 
writeln(Output,Mission[15] 
writeln(Output,Mission[16] 
writeln (Output ,Mission[17] 
writeln (Output) ; 


:10,/CE-102G 2 

2:10, LeadTime:10,INSTL:10) ; 
210, ST: 10, Hor ue. 
:10,Meth:10) ; 

:10,Rep:10) ; 


end; 
{+++PART2 ESTIMATE THE DEGRADATION+++++++++++++++++++++++} 
{ These procedures estimate the time of onset of subsystem 


degradation and the magnitude 
over time. It includes: 
ReadData 
BuildVec 
Mu_Eta 
SigmaSqr_Sum 
FindHat Parameterl 
InitialSigmaSqr 
FindHat Parameter2 


Procedure ReadData 
var XVec 
var YVec 

var I integer; 

F, MA, TH real; 

{ This procedure reads data 

time between failure 


(var Input 


(maintenance action) 


and evolution of the degradation 


(both method), 

(both method), 

(maximum likelihood), 

(Bayesian), 

(Bayesian). } 


text; 
RealVec; 
RealVec) ; 


from a file, and transforms mean 
to number of 


failures (maintenance actions). } 
Procedure FindLength (Var Input text; 
var I integer) ; 
var F, MA real; 
begin 
reset (Input); 
i =r 
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while not eof (Input) do 


begin 
while not eoln (Input) do 
begin 
read(Input,F,MA); 
tos = t+. 
end; 
readln (Input); 
end; 
Pees = I; 
end; 


LHF ttt tt Ht tte tt tettetttettttete+ete+tt+ i} 


Procedure Compare (XVec, YVec : RealVec; 


I : integer) ; 
begin 
if XVec[I] >= YVec[I] then 
begin 
if XVec[I]) > Max then 
Max := XVec[I]; 
if YVec[I] < Min then 
Mon := Yweelr]; 
end 
else 
if YVec[I]) > Max then 
Max := YVec[I]; 
if XVec[I]) < Min then 
Min := XVec[I]; 
end; 


LHEEEHEEEHE HEHEHE HE H+ ttt tttttet+ttt++t++tts i} 


begin {Procedure ReadData} 
assign (Input,Infile) ; 
FindLength(Input,I); 
reset (Input) ; 
hie s= NS~UPMe 
while not eof (Input) do 
begin 
while not eoln (Input) do 
begin 
read (Input,F,MA); 
if F <> 0 then 


KVec(I) := TH/F 
else 

KXVec{I] := Os 
if MA <> 0 then 

YVec[I] := TH/MA 
else 
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Yvec (i= 0; 
Compare (XVec, YVec,I); 
l := t=, 
end; 
readln(Input) ; 
end; 
close (Input); 
end; {Procedure ReadData} 


LEEEEEEEEELEEEE EEE EERE EE EE EEE EEE EEE EEE FEE HEF t+ t+ tt tettttet+++t++++ i 
{ These procedures use the given parameters (SetMuF, 
SetEtaF, SetSigmaSqrF, SetCF, SetMuMA, SetEtaMA, 
SetSigmaSqrMA, SetCMA) to generate two simulating data set. 
All parameters are used to generate number of failures 
(maintenance actions) instead of mean time between failure 
(maintenance action). } 


Function GenNormal : real; 
var Uni : real; 


Function GenUniform : real; 
var U: real; 


begin {Function GenUniform} 


U := Random; 

if U > 0 then 
GenUniform := U 

else 
GenUniform := Random; 


end; {Function GenUniform} 


begin {Function GenNormal} 

Uni := GenUniform; 

GenNormal := Sqrt((-Ln(Un1i) *2) )*Cos(2*Pi*Random) ; 
end; {Function GenNormal} 


LEEEEEEFHEE EE EE EEF EE EE FETHEH H+ +++ ++ +444} 


Procedure BuildVec (T, SetC : integer; 
SetSigmaSgqr, SetMu, SetEta : real; 
var XVec : RealVec) ; 
var I, DF, D: integer; 
NorF, Nor, SigmaF, Sigma : real; 


begin {Procedure BuildVec} 


Sigma := sqrt(SetSigmaSqr) ; 
for L, -= eto, Teldec 
begin 
Nor := GenNormal; 
if I <= SetcC then 
Da = 0 
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else 
Deel =Secc- 
XVec[I] := SetMu+(D*SetEta)+(Sigma*Nor) ; 
end; 
end; {Procedure BuildVec} 


{+++ +++ 44+ 4444+ +444 4444444444444 4+4+4+44+4+4+} 


Procedure BuildVecl (T, SetC : integer; 


SetSigmaSqr, SetMu, SetEta : real; 
var XVec : RealVec); 
var I, DF, D: integer; 
NorF, Nor, SigmaF, Sigma : real; 


begin {Procedure BuildVec1} 
Sigma := sqrt (SetSigmaSqr) ; 
Meet <= 1 to T de 
begin 
Nor := GenNormal; 
if I <= SetcC then 
pa s= © 
else 
DD -=si-Seee- 
XVec[I] := SetMu+(D*SetEta)+(Sigma*Nor) ; 
if XVec[I] > Max then 
Max := XVec[T] 
else 
if XVec[I]) < Min then 
Min := XVec(I]; 
end; 
end; {Procedure BuildVec1} 


LHEFEEEHE HEHEHE FE FHE EH EHEE FEE FHF HE ++ 4+ +++ tetettettetetttetet+t+etteet4)} 
{ This procedure calculates the mean value of array before 


the trend, the occuring time of the trend, and the slope of 
the trend. } 


Procedure Mu_Eta (T, C : integer; 
var Phil : real; 
var Phi2 : real; 
XVec : RealVec; 
var Mu : RealVec; 
var Eta : RealVec); 

var Xl, X2, Phi: real; 


{HEHEHE HEHEHE HEHE HH +H teetteeteett+et++t+ } 
Function BarXl (T : integer; 
XVec : RealVec) : real; 


var I > integer; 
Sum : real; 
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begin {Function Barx1} 


Sum :=30F 

for I := 1 tome 
Sum := Sum+XVec[I]; 

Barxl := Sum/7@- 


end; {Function BarXx1} 
LHFEEEEHHEE HEHEHE HH HHH ++ tee teetett++t+t+} 
Function BarxX2 (T, C : integer; 
XVec : RealVec) : real; 
Var I : integer; 
Sum : real; 


begin {Function BarxX2} 


Sum := 9; 
for 0 s= 7 (C+1) to, Laae 

Sum := Sum+(XVec[I]*(I-C)); 
BarxX2 := Sum/T; 


end; {Function Barx2} 

LHFHEHHEHHEHHEHEAHHE TH+ H+ + +++ te ttet++t+t++++++i 
Function BarPhil (T, C : integer) : real; 
Var 1 : integer; 


Sum : real; 


begin {Function BarPhil} 


Stm = 0; 

for Ees=—1 tosh eimds 
Sum := Sum+I; 

BarPhil := Sum/T; 


end; {Function BarPhil} 

{HEHEHE HE HEHEHE +H HHH +H +++ t+eetetet+et+t++t+++++) 
Function BarPhi2 (T, C : integer) : real; 
var I : integer; 


Sum, Temp : real; 


begin {Function BarPhi2} 


Sim 

for I <= 1 tomer -C)} de 
Sum := Sum+(I*I); 

BarPhi2 := Sum/T; 


end; {Function BarPhi2} 
LHEEEEH HEHEHE EHEEEE HEH HEEE HE +H+ HE +++ ++¢++4+4+4+4} 
begin {Procedure Mu_Eta} 


X1 := Barxl(T,xXVec); 
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tae < | then 


begin 
X2 = Barx2 (T,C, XVec) ; 
Phil = Barend ec); 
Phi2 = BarPhi2(T,C); 
Phi = Phi2-(Phil*Phil); 
Mu[C] = ((Phi2*X1)-(Phil1*X2))/Phi; 
Eta[C] := (X2-(Phil*Mu([C]))/Phi2; 
end 
else 
begin 
Male) 2= Xl: 
Beale := 0; 
end; 
end; {Procedure Mu_Eta} 


LEA FEEEEEEFEPEEEPEEEEEE EEE EEE EEE HEHEHE F EEE HEH +++ + 44444444444} 


Function SigmaSqr_Sum (T, C : integer; 
XVec, Mu, Eta : RealVec) : real; 
var I, D : integer; 
Temp, Sum : real; 


begin {Function SigmaSqr_Sum} 


Sum := 0; 
bor! -= 1 towel de 
begin 
if I <= C then 
pee = 0 
else 
Dee = 1 -C); 
Temp := XVec[1I]-Mu[(C)-(D*Eta[C] ); 
Sum := Sum+(Temp*Temp) ; 
end; 
SigmaSgqr_Sum := Sum; {Likelihood : /T, Bayesian 
yo 1 } 


end; {Function SigmaSqr_Sum} 
{FEAF EEEEE EEE FE EE EEEEEE EE EE EERE EEE EEEEE EEE EH ETH tHe e+ 4st } 


Procedure FindHatParameterl(T : integer; 
var HatC : integer; 
XVec : RealVec; 
var Mu : RealVec; 
var Eta : RealVec; 
var SigmaSqr : RealVec; 
var HatMu : real; 
var HatEta : real); 

var C :integer; 

Phil, Phi2, S, MinS, HSS : real; 


o1. 


begin {Procedure FindHat Parameterl1 } 
MinS := 0; 
for C := 1 femigaeo 
begin 
Mu_Eta(T,C, Phil, Phi2,XVec,Mu, Eta) ; 
SigmaSqr[C] := SigmaSqr_Sum(T,C,XVec, Mu, Eta) /T; 
1£ SigmaSqr[C] > 0 then 
S := l+In(Si10qmacar lew 
else 
S := -1E+10; 
i1f (MinS = 0) or (MinS >= S) then 
begin 
MinS 
Hat¢€ 
end; 
end; 
HatMu := Mu[HatC]; 
HatEta := Eta[HatC]; 
end; {Procedure FindHatParameter1} 


oF 
a 


LEAF EEEEHEEEEFEFEEEFHEEE HEE EEEFEEFEFEFE FEE HEHEHE E EEE + 444444+4+4444} 


Procedure InitialSigmaSqr (XVec : RealVec; 
var Mu : RealVec; 
var Eta : RealVec; 
var SigmaSqr : RealVec); 
var I : integer; 
Sum : real; 


begin {Procedure InitialSigmaSqr} 

Sum s= 0; 

for I := 1 toes de 

Sum := Sum+XVec[I]; 

Mal5] >= sSumy 5; 

Pears] :="0; 

SigmaSqr[5] := SigmaSqr_Sum(5,5, XVec,Mu, Eta) /4; 
end; {Procedure InitialSigmaSqr} 


LEFF EEEEAEEEEEEE FE EEE EE EEE FEFEFEAFA FE FHEEEHE EEE ++++44+44+44+4} 


Procedure FindHatParameter2 (T : integer; 
XVec : RealVec; 
var Mu : RealVec; 
var Eta : RealVec; 
var SigmaSqr : RealVec; 
var VSqr : RealVec; 
var RSqr : RealVec; 
var RhoSqr : RealVec; 
var PiStar : RealVec) ; 

var C : integer; 

Sum, Suml, Sum2 : real; 
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Phil, Phi2, Value : real; 
Coef, Sig2, KStar : RealVec; 


LEAH EEEEEEEE EEE EEE EHH HHH +44 +tetteett} 


Procedure V_R_RhoSqr (T, C : integer; 
Phil, Phi2 : real; 
SigmaSqr : RealVec; 
var VSqr : RealVec; 
var RSqr : RealVec; 
var RhoSqr : RealVec; 
var Coef : RealVec); 

var Phi : real; 


begin {Procedure V_R_RhoSqr} 
me C < T then 


begin 
Phi := Phi2-(Phil*Phil); 
VSqr[C} = SigmaSqr[T-1]/T/Phi; 
RSqr [C] = VSar ie] *Pha2: 
RhoSqr[C] := Phil*Phil/Phi2; 
Coef[C] =e biccante (Cl Rhos |e) je Rsaqr(C)]*vSqr[c)); 

end 

else 

begin 
RSqr[C] = SigmaSqr[T-1]/T; 
Coef[C]} = Samet (2*baeRSqr([C]); 
Veqr[(c] = 0; 
Fhocar(C] := 0; 

end; 


end; {Procedure V_R_RhoSqr} 
LEER EEHEEEEE HEHEHE EFF HEHEHE +++ ++tt+++++4+4+4 i} 


Function Beta (K, T : real): real; 
var N, D: real; 


{ This function computes the prior function combined with 
geometric probability function. 

It is given that both coefficients of beta (prior) 
function are 1. } 


begin {Function Beta} 
if K = 1 then 
Beta := 0.5 
else 
if K =2 then 
Beta := 1/6 
else 
mroke<= T then 
Beta := 1/(K* (K+1)) 
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else 
Betas) K> 
end; {Function Beta} 


LHFEFEHHEAHEE FE HFHEFEFF EFF +4 4+44+4+4¢4+4+4+4+4+4+4+4+4} 


begin {Procedure FindHatParameter2} 
Sun ee 
for C := 1 to (iri ede 
begin 
Mu_Eta(T,C, Phil, Phi2, XVec atu, Etsa) ; 


V_R_RhoSqr (T,C, Phil, Phi2,SigmaSqr,Vsqr, RSqr, Ra@eenr, Gooam 
Sig2Z[€] := SigmaSqr_Sum(@T, CC, .Vec Mur eea 
KStar(C] := Sig2[C]/(2*SigmaSar(T-1]); 

Value := 50-KStar([C]; 
if Value >= -86 then 

PiStar[C] := exp(Value) *Beta(C, Tf) *Ceee [C] 
else 

PiStar[C]) := 0; 

{ PiStar[C] := Beta(C,T) *exp(-KStar[C] ) *Coef[C];} 
Sum := Sum+PiStar[C]; 

end; 


Suml := Q; 
forve 2= loco (lee) aae 
begin 
PiStar[C] := PiStar [€j}]/sum; 
if PiStar[C]) < 1E-6 then 
Pistar[C] := 0; 
Suml := Suml+PiStar([C]; 
end; 


Sum2 := 0; 
tor €C <= 1 to (eeligde 
begin 
PiStar[(C] := PiStar[C])/Suml; 
Sum2 := Sum2+(PiStaxy[C] *Saq2iie)/-f) ie 
end; 
SigmaSer[T] 7="Sum2; 
end; {Procedure FindHatParameter2} 


{+++PART3 : COMPUTE THE COSTS & SEARCH BEST CHANGE TIME ++++} 
{ These procedures compare the cost of remaining with the 
current subsystem with the cost of investing in the upgraded 
subsystem to obtain a best time to invest in the upgraded 
subsystem. It includes: 


HatCost_New (both method), 
HatCost_NewAD (both method), 
HatCost_Old (both method), 


HatCost_OldAD (both method), 
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FindCostNew (Bayesian), 


FindCost0Old (Bayesian), 
Mean_Variance (Bayesian), 
SearchTime (both method), 


FindChangeTimel (maximum likelihood), 
FindChangeTime2 (Bayesian), 
Uncertaintyl (maximum likelihood), 
Uncertainty2 (Bayesian). 


Procedure Install (H, T, Tau : integer; 
var N : integer; 
var PartialSum : real); 


begin {Procedure Install} 
Pe s= trunc (NS/INSTB) ; 
if N > (H-(T+Tau)) then 
N := H-(T+Tau) ; 
PartialSum := (N* (N+1)/2)*(INSTL/NS) ; 
end; {Procedure Install} 


LEFF EEEEEEEEEEEEEEEEHEEFEEF+E++4444+4+4++4} 
Function HatCost_NewAD (H, T, Tau, N : integer; 
PartialSum, CO, CN : real): real; 


var NewPart, OldPart : real; 


begin {Function HatCost_NewAD} 


NewPart := CN*(PartialSum-N+H- (T+Tau) ); 
OldPart := CO* (Tau+1+N-PartialSum) ; 
HatCost_NewAD := OldPart+NewPart; 


end; {Function HatCost_NewAD} 

{FEE EEEHE HEHEHE FEE HEHEHE HEHEHE ++ 44+4+4+4+4+4+4+4+4+4+4} 

Function HatCost_New (H, T, Tau, HatC, N : integer; 
PartialSum, HatMu, HatEta : real; 
CO, CN, NONew : real): real; 


var S, D: integer; 
Sum, Partial, NewPart, OldPart : real; 


begin {Function HatCost_New} 


NewPart := (CN*NONew) * (PartialSum-N+H- (T+Tau) ); 
if HatC >= T then 

OldPart := (CO*HatMu) * (Tau+1+N-PartialSum) 
else 
begin 

® := T-HaecC; 

Sum <= 0% 

for S -= 0 to (Tau+N)do 


25 


i1£°98 <= Tau then 


Partial := 1 
else 
Partial := 1-(S-Tau) * (INSTL/NS) ; 
Sum := Sum+(HatMu+HatEta* (S+D))*Partial; 
end; 
OldPart := CO*Sum; 
end; 


HatCost_New := OldPart+NewPart; 
end; {Function HatCost_New} 


LH FHEEHEE EEE EH FEEFE EEE H HEHEHE HHH ++ t+ tet+tets i} 


Function HatCost_OldAD (H, T : integer; 
CO : real) : real; 


begin {Function HatCost_Ol1dAD} 
HatCost_OldAD := CO*(H-(T-1) ) 
end; {Function HatCost_Ol1dAD} 


LEAF EEEEEE HEHEHE HEHEEFEFHEFE +++ + +4444+4+4+44+4} 


Function HatCost_Old (H, T, HatC: integer; 
HatMu, HatEta, CO : real) : real; 
var S, D : integer; 
Sum : real; 


begin {Function HatCost_0Old} 
if HatC > T then 
HatCost_Old := CO*HatMu* (H-(T-1) ) 


else 
begin 
D:= TSHace: 
Silja 0: 
for S := 0 to (H=wiedo 
Sum := Sum+HatMu+HatEta* (S+D) ; 
end; 


HatCost_Old@ <;= CO™sum- 
end; {Function HatCost_0Ol1d} 


LEFF FEE EFEFE FE FEAF EE EEEEAFEFEE HEE EE FEA EEF E EFF +A +++ 44444444444} 


Function FindCostNew (H, T, Tau, N : integer; 
PartialSum, CO, CN, NONew : real; 
Mu, Eta : RealVec; 
PiStar : RealVec) : real; 
var Sum, HCN, HatMu, HatEta : real; 
C : integer; 


begin {Function FindCostNew} 
Sune := 07 


a6 


fer C := 1 to (T-LbeadTime+1) do 


begin 
HatMu := Mu([C]; 
HatEta := Eta[C]; 
af Pastas (Cle<> 0 then 
HON := HatCost_New (H, T, Tau, C)N; PartialSum, 
HatMu, HatEta,CO, CN, NONew) *PiStar([C] 
else 
ery := 0; 
Sum := SumrHCn; 
end; 
FindCostNew := Sum; 


end; {Function FindCostNew} 
{HEHEHE HEHEHE HHH +t tee ett tee et t44+t444++ i 


Function FindCostOld (H, T : integer; CO : real; 
Mu, Eta : RealVec; 
PiStar : RealVec) : real; 
var Sum, HCO, HatMu, HatEta : real; 
C : anteger; 


begin {Function FindCost0ld} 


Sum := 0; 
for C := 1 to (T-LeadTime+1) do 
begin 

HatMu := Mu[C]; 

Haebea :=eBta(C]> 


TF" PiStar(c]) <> 0 then 
HCO := HatCost_Old(H,T,C,HatMu, HatEta,CO) *PiStar([C] 


else 
HEO :=ea0; 
Sum := Sum+HCoO; 
end; 
Fin@cesteld := Sum-: 


end; {Function FindCost0Old} 
LEER HEHE HEEHEE EE EEE EEE EE EEE HEHEHE EEE HEHEHE HEHEHE HEH ee et tt teste tte+t+s+)} 


Procedure Mean_Variance (H, T, Tau, N : integer; 
PartialSum, CO, CN, NONew : real; 
Mu, Eta : RealVec; 
Roar, Vaan, Rhosaqr : RealVec; 
PiStar : RealVec; 
var CostMean : real; 
var CostWaxr2ance : real); 
var C : integer; 
Sum, CostVanl, CostVarn2 : realy; 
MuCo, ConstCo : real; 
EtaCo, CondMean : RealVec; 


oF 


{HAHAH HEHEFHHE +++ +++ 44+44+444+44+4+4+4+4+4+44} 


Function ConstCoef (H, T, Tau, N : integer; 
PartialSum, CN, NONew : real): real; 


begin {Function ConstCoef} 
ConstCoef := (CN*NONew) * (PartialSum-N+H- (T+Tau) ); 
end; {Function ConstCoef} 


{+++ H+++HH+E HHH H +++ +44+4+44+444+4+4+444+44+444 } 
Function MuCoef (H, T, Tau, N : integer; 
CO : real): real; 
var S : integer; 
Sum, Partial : real; 


begin {Function MuCoef} 


Sunes =a 
for S := 0 tov(TausN) de 
begin 
if S <= Tau then 
Partial := 1 
else , 
Partial := 1-(S-Tau) * (INSTL/NS) ; 
Sum := Sum-Partial; 
end; 
for S := 0 Een (h-)) deo 
Sum := Sum+l1; 
MuCoet := CO*Stm; 


end; {Function MuCoef} 
{HAEEEEHHEHEE HEHEHE HHP HE +++ +4+44+4+4+4+4+4+4+4+4} 


Function EtaCoef (H, T, Tau, HatC, N : integer; 
CO : real): real; 


var S, D: integer; 
Sum, Partial : real; 


begin {Function EtaCoef} 
if HatC >= T then 


EtaCoef := 0 
else 
begin 
D := T-HatGy 
Sum. 7=e0- 
for S := 0 to (Tawen) do 
begin 
if S <= Tau then 
Partial := 1 
else 


a8 


Partial := 1-(S-Tau) * (INSTL/NS) ; 


Sum := Sum-(S+D)*Partial; 
end; 
Eoceom:= 0 to (H-T) do 
Sum := Sum+(S+D); 
end; 
Bitacoef := CO*Sum; 


end; {Function EtaCoef } 
LHF HHEHEH HEH HHH HHH + tee tee te+tetts i} 


begin {Function Mean_Variance} 

Muco := MuCoef (H,T,Tau,N,CO) ; 

ConstCo := ConstCoef(H,T,Tau,N, PartialSum, CN, NONew) ; 

sum := O> 

for C := 1 to (T-LeadTime+1) do 

1f PiStar[C] <> 0 then 

begin 
EtaCo[C] := EtaCoef (H,T,Tau,C,N,CO); 
CondMean[C] := MuCo*Mu[C]+EtaCo[C]*Eta[C]-ConstCo; 
Sum := Sum+CondMean[C]*PiStar[C]; 

end; 

GeseMean := Sum; 


Sum := 0; 
for C := 1 to (T-LeadTime+1) do 
1£ PiStar[C] <> 0 then 
Sum := Sum+sqr(CondMean[C]-CostMean) *PiStar[C]; 
CesceVarZ == Sum; 


Sum 2= 0; 
for C := 1 to (T-LeadTime+1) do 
1f PiStar[C] <> 0 then 
Sum := Sum+(sqr(MuCo) *RSqr[C]+sqr(EtaCo[C]) *VSqr[C] 
=2* Sqr ERNOSar [(@} *RSaqrifie] *VSqr [(C]}9 
*MuCo*EtaCo[C])*PiStar([C]; 
CostvVaerl :=e«Sum; 


CostVariance := CostVarl+CostVar2; 
end; {Function Mean_Variance} 


LEER EEE EERE EEE FE EEE FEE FEE EEE FEE EEE EEE EEE FEE FETE +++ 444444} 


Procedure SearchTime (H, N : integer; 
HC : real; 
var Tau : integer; 
var CT : integer; 
var MinC : real); 


begin {Procedure SearchTime} 
if HC > MinC then 
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begin 
if MincC > 0 then 


CT -=eh 
else 
CT 2= Gila: 
Tats: =sen 
end 
else 
begin 
MincC := HC; 
1f CT = (H-LeadTime-N) then 
begin 
if Minc > 0 then 
Gro: =H: 
Taw v= H; 
end; 
end; 
end; {Procedure SearchTime} 


LFEEEEEEEEFEFEEEEFEEEEEEEEFEEEEEEEE EEE HEHEHE +++ +++ 4+44+4+4+4+4+4+4+4++4} 
{ Procedures "FindChangeTimel" and "FindChangeTime2" 
estimate the best time for subsystem upgrade. } 


Procedure FindChangeTimel (H, T, HatCF, HatCMA : integer; 
HatMuF, HatEtaF : real; 
HatMuMA, HatEtaMA : real); 
var N, Tau, CT, Tl : integer; 
HC, MinC, PartialSum, CostNew, CostOld, CostFix : real; 


begin {Procedure FindChangeTimel } 


Tait := X). 

Tl := T-LeadTime; 

Install (H,T, Tau,N, PartialSum) ; 

CostFix := CF; 

CostNew := HatCost_New(H,T,Tau,HatCF,N, PartialSum, 


HatMuF, HatEtaF, COF, CNF, NOFNew) 
+HatCost_New(H,T, Tau, HatCMA,N, PartialSum, 

HatMuMA, HatEtaMA, COMA, CNMA, NOMANew) 
+HatCost_NewAD (H,T, Tau,N, PartialSum, COAD, CNAD) 
+COStFix; 

CostOld := HatCost_Old(H,T,HatCF, HatMuF, HatEtaF,COF) 
+HatCost_Old(H,T,HatCMA, HatMuMA, HatEtaMA, COMA) 
+HatCost_OldAD(H,T,COAD) ; 

MincC := CostNew-CostOld; 


if MinC > CF then 
CTvVec(T1] sen 

else 

begin 
while Tau <= (H-T) do 
begin 
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Tau :=8Tau+1- 
eltee=) Tle Tay: 
if CT <= (H-LeadTime-N) then 
begin 
if H = (T+Tau) then 
CostFix := 0; 
CostNew := HatCost_New(H,T, Tau, HatCF,N, PartialSum, 
HatMuF, HatEtaF, COF, CNF, NOFNew) 
+HatCost_New(H,T,Tau,HatCMA,N, 
PartialSum, HatMuMA, HatEtaMA, 
COMA, CNMA, NOMANew) 
+HatCost_NewAD (H,T,Tau,N, PartialSum, 


COAD, CNAD) 
+CostFix; 
HC := CostNew-CostOld; 
SearchTime(H,N,HC,Tau,CT,MincC) ; 
end 
else 
ee —o8) 5 
end; 
CTVeewiT1] := CTy 


end; 

write(T1:5,CTVec(T1]:8); 

mierte (Output, T1:5,CTVec([T1]:8) ; 
end; {Procedure FindChangeTimel } 


LEFEEEEEEEFEEEEEEEFEEEEEEE FEET FEET EEE EEE FE EEE FEE H+ + ++4+4+4+4+4} 


Procedure FindChangeTime2 (H, T : integer; 
MuF, EtaF : RealVec; 
MuMA, EtaMA : RealVec; 
PiStarF, PiStarMA : RealVec); 
var N, Tau, CT, Tl : integer; 
HC, MincC, PartialSum, CostNew, CostOld, CostFix : real; 


begin {Procedure FindChangeTime?2 } 


Paes = 0: 

Tl := T-LeadTime; 

Install (H,T,Tau,N, PartialSum) ; 

CostFix := CF; 

CostNew := FindCostNew(H,T, Tau,N, PartialSum, 

COF, CNF, NOFNew, MuF, EtaF, PiStarF) 
+FindCostNew(H,T,Tau,N, PartialSum, 

COMA, CNMA, NOMANew, MuMA, EtaMA, PiStarMA) 
+HatCost_NewAD (H,T, Tau,N, PartialSum, COAD, CNAD) 
+CostFix; 

CostOld := FindCost0Old(H,T,COF,MuF,EtaF, PiStarF) 
+FindCostOld(H,T,COMA,MuMA, EtaMA, PiStarMA) 
+HatCost_OldAD(H,T,COAD) ; 

MincC := CostNew-CostOld; 


AOE 


if MinC > CF Ghen 


CTVec(Ti 
else 
begin 
while Tau <= (H-T) do 
begin 
Tau := PTauyl: 
CT <-= Ti-Tau- 
if CT <= (H-LeadTime-N) then 
begin 
if H = (T+Tau) then 
CostFix := 0; 
CostNew := FindCostNew(H,T,Tau,N, PartialSum,CoOF, 
CNF, NOFNew, MuF, EtaF, PiStarF) 
+FindCostNew (H,T, Tau,N, PartialSum, COMA, 
CNMA, NOMANew, MuMA, EtaMA, 
PiStarMA) 
+HatCost_NewAD(H,T,Tau,N, PartialSum, 
COAD, CNAD) 
+CostFix; 
HC := CostNew-CostOld; 
SearchTime(H,N,HC,Tau,CT,Minc) ; 
end ‘ 
else 
CT. She 
end; 


efvVec?rl)—:="Cr; 
end; 

write(T1:5,CTVec[T1] :8); 

write (Outpubey TIS, CTVec (Thm). 
end; {Procedure FindChangeTime2 } 


LHAEEE EEA EE FEEEE EEE EEE EEE EEE EE EE EE EE EEE EEE THAT Ht +++ tt++++t++) 


{ Procedures “Uncertaintyl" and "Uncertainty2" compute the 
MEAN, SD, and 2 SD bounds of cost advantage of upgrade for 
the estimated best upgrade time. } 


Procedure Uncertaintyl (H, T, R, HatCF, HatCMA : integer; 
HatMuF, HatEtaF : real; 
HatMuMA, HatEtaMA : real); 
var PartialSum, CostNew, CostOld, Sum, SD : real; 
N, Tau : integer; 


begin {Procedure Uncertaintyl} 


Tau := CTVec[{Time] -Time; 
if CTVec{Time] = Hor then 
Tau := 0; 


CostOld := HatCost_Old(H,T,HatCF, HatMuF, HatEtaF, COF) 
+HatCost_Old(H,T,HatCMA, HatMuMA, HatEtaMA, COMA) 
+HatCost_OldAD(H,T,COAD) ; 

Install (H,T,Tau,N, PartialSum) ; 
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CostNew := HatCost_New(H,T,Tau, HatCF,N, PartialSum, 
HatMuF, HatEtaF, COF,CNF, NOFNew) 
+HatCost_New(H,T, Tau, HatCMA,N, PartialSum, 
HatMuMA, HatEtaMA, COMA, CNMA, NOMANew) 
+HatCost_NewAD(H,T,Tau,N, PartialSum, COAD, CNAD) 
+CF; 
Gain[R] := CostOld-CostNew; 
Total := Total+Gain[R] ; 


if R = Rep then 
begin 
sum += 0% 
Mean[Time] := Total/Rep; 
forgeee= 1 “bor ReEp do 
Sum := Sum+sqr(Gain[R]-Mean[Time] ); 
SD := sqrt(Sum/(Rep-1)); 
LB[Time] := Mean[Time] -2*SD; 
UB[Time] := Mean[Time]+2*SD; 
tfetin > LB[Time) then 
Min := LB[Time]; 
if Max < UB[Time] then 
Max := UB[Time]; 


writeln(’ ‘,Mean[Time]:10,’ ‘, 
SBei0),.’ ', UB [tame |S 1:09’ ‘, UB[Time] :10); 
writeln(Output, ’ ’ MeamiTime):10,° ‘, 
SpD:10,' *’, LB([Timey: 10,’ ‘,UB([Time] :10); 
end; 


end; {Procedure Uncertaintyl} 
LFEEEEEEEH EEE EE EE EEEE EEE EEE EEE HEHEHE HEF HEHEHE THE HH +H + t4++4+44+4+4+4} 


Procedure Uncertainty2 (H, T : integer; 
MuF, EtaF, MuMA, EtaMA : RealVec; 
RSqrF, VSqrF, RhoSqrF : RealVec; 
RSqrMA, VSqrMA, RhoSqrMA : RealVec; 
PiStarF, PiStarMA : RealVec) ; 
var PartialSum, MeanF, MeanMA : real; 
VarianceF, VarianceMA, SD : real; 
N, Tau : integer; 


begin {Procedure Uncertainty2} 


Tau := CTVec[Time] -Time; 
if CTVec[Time] = Hor then 
baie: = > 


Install (H,T,Tau,N, PartialSum) ; 

Mean_Variance(H,T,Tau,N, PartialSum, COF, CNF,NOFNew, 
MuF,EtaF,RSqrF,VSqrF,RhoSqrF, PiStarF, 
MeanF, VarianceF) ; 

Mean_Variance(H,T,Tau,N, PartialSum, COMA, CNMA, NOMANew, 
MuMA, EtaMA, RSqrMA, VSqrMA, RhoSqrMA, PiStarMA, 
MeanMA, VarianceMA) ; 
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SD := Sqrt (VarianceF+VarianceMA) ; 

Mean[Time] := MeanF+MeanMA-CF 
-HatCost_NewAD (H,T, Tau,N, PartialSum, COAD, CNAD) 
+HatCost_Ol1dAD(H,T,COAD) ; 


LB[Time] := Mean[Time] -2*SD; 
UB[Time] := Mean[Time]+2*SD; 
if Min > LB[Time] then 

Min := LB[Time] ; 


if Max < UB[Time] then 
Max := UB[Time] ; 


writeln(’ ‘,Mean[(Time]:10,' ‘, 
SD:10,’ ‘,LB[Time):10,' ‘,USB(Time] eto 
writeln(Output, ’ ‘,Mean[Time):10,’ ’, 
SD:10,' ‘,LB[Time]:20,° Svusiainel 
end; {Procedure Uncertainty2} 


LAF EEEEEEEEEFEAEEEEEF EEE EE FEET EE EEE EEE EFEFHEE+E HEH 4+ 444444444} 
{ These procedures reset the values for assessing 
uncertainty. 

1 for maximum likelihood, 

2 for Bayesian. } 


Procedure ResetParameterl1; 5 


begin {Procedure Reset Parameter1} 
SeCtCF > =e8Haeer- 


SetMuF := HatMuF; 
SetEtaF := HatEtaF; 
SetSigmaSqrF := SigmaSqrF[SetCF] ; 


SetCMA := HatCMA; 

SetMuMA := HatMuMA; 

SetEtaMA := HatEtaMA; 

SetSigmaSqrMA := SigmaSqrMA[SetCMA] ; 
end; {Procedure ResetParameter1} 


LFEEEEALEFEFE HEHEHE HEE FEHHEFEF+EF +++ +44+4+4+4+4+4} 


Procedure ResetParameter2; 
var T : integer; 


begin {Procedure Reset Parameter? } 


PiDistF[1] := PiStarF[1]; 
for T := 2 to (Time+l1) do 
PiDistF[T]) := PiDistF[T-1)+PiStarF[T] ; 
PiDistMA[1] := PiStarMA[1]; 
for T := 2 to (Time+l) Go 
PiDistMA[T] := PiDistMA[T-1]+PiStarMA[T] ; 
end; {Procedure Reset Parameter2 } 
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{+++PART4 : GENERATE DATA GRAPH 4++++4++4++++++++++++++++++++++++} 
{ These procedures plot the graph of the mean number of 


failures and mean number of maintenance actions. 


Procedure GenDataGraph; 
var GraphDriver, GraphMode : integer; 
YScale, XScale : real; 


{HEE HEHEHE HEHEHE HEHEHE HHH ete tttetetteteett+++ } 


Procedure DrawAxis; 
var S : String; 


begin {Procedure DrawAxis} 
SetColor (White) ; 
Out vextxy(10,10,' NO. of Failures’); 
OutTextXY(10,20, ‘(Maintenance Actions) ’); 
Out Text XY (250,10, ‘Mean number of Failures’); 


Out Text XY (250,20, ‘Mean number of Maintenance Actions’); 


faneis0,30,30,430); 
mone (30,430,630,430); 
end; {Procedure DrawAxis} 


LHEEEEEEEEEHEEEEE HEHE FE HEE HHH +4 44444444444} 
Procedure SetYScale (YScale : real); 
var I, Y : integer; 

M, Scale, Length, Units : real; 

See string; 


begin {Procedure SetYyScale} 


he 2= 10> 
while (Max/M) > 10 do 
M := M*10; 
Scale := M/2; 
Length := 0; 
Units := 0; 
while Length > -Min do 
begin 
Length := Length-Scale; 
Units := Units+0.5; 
1f"Units > VO then 
units := Units-10; 
end; 
I := round(frac(Units)) ; 
while Length > -Max do 
begin 
Y := round((Length+Min) *YScale) ; 
Line (30, 430+Y, 630, 430+Y) ; 
wc. (bemod 2) = 0 then 
begin 
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Str (Unzess i 0esie 
Out Text XY (15, 430+Y,S); 


end; 
Length := Length-Scale; 
Units := Units+0.5; 
te I+1; 

end; 


Str(round (hs) 
Out Text XY (30,470, 'Y-Axis scale’); 
Out Text Xy (130749707 S) ; 

end; {Procedure SetYScale} 


LHAFEEHEHEHE HEHEHE EHH H+ +++ +++ ttte+t+++4+4+4} 
Procedure SetXScale (XScale : real); 
var T, Tl, X : integer; 


Sse rng, 


begin {Procedure SetXScale} 


== 5 
while T <= ET do 
begin 
X := round((T-1)*xScaley 


Line (30+X,30,30+X, 430); 
if (T mod 10) = 0 then 
begin 
Str (13s) 
Out Text XY (15+X,440,S); 
end; 
=e eo 
end; 
Out Text XY (300, 460, ’Month’ ); 
end; {Procedure SetXScale} 


LA FEEEE HEHEHE EEE HEEHEEFEA HEHEHE +A ++ +4 4444444 } 


Procedure DrawLine (YScale, XScale, Mu, Eta : real; 
C : integer); 
var Xi, Y1, 92 <= integer; 


begin {Procedure DrawCurve} 
X1 := 30+round((C-1)*XScale); 
Y¥1 := 430-round( (Mu-Min) *YScale) ; 
bine (30,71. ae Vie 
if ET > € then 
begin 
Y2 := 430-round( (Mu+Eta* (ET-C) -Min) *YScale) ; 
Line (X11, Y1, 6307 Y20e 
end; 
end; {Procedure DrawCurve} 
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{++ +H +++ +++ +++ ++ +++ ++ ++t+tt+tsettet+t++s i} 


Procedure DrawCurve (YScale, XScale : real; XVec : RealVec); 
var Xl, X2, T : integer; 
YM1, YM2 : integer; 


begin {Procedure DrawCurve} 
ce nl 
el aataae 5)() : 
YM1 := 430-round( (XVec[T] -Min) *YScale) ; 
moe] !=— tO ET do 
begin 
X2 := round( (T-1) *Xscale+30) ; 
YM2 := 430-round((XVec[T] -~Min) *YScale) ; 
Line (X1,YM1,X2, YM2) ; 
Gui rextxy (X2-3, YMZ2-5,'’.'): 
me = K2> 
WM 8 := YM2; 
end; 
end; {Procedure DrawCurve} 


{+4+4+44+44+444+4¢4+4¢+4+44+4+44++4+4+4+4+4+4+4+4+44+44} 


begin {Procedure GenDataGraph} {Verify Driver, Mode, and path} 
GraphDriver := VGA; {adjust to the acceptable values} 
GraphMode := VGAHi; 
InitGraph(GraphDriver, GraphMode, ‘'a:\bgi\'); 
Set BKColor (Black) ; 


ClearDevice; 
DrawAxis; 
YScale := 400/(Max-Min) ; 


XScale := 600/(ET-1); 
SetLineStyle(DottedLn, 0,NormWidth) ; 
SetYyScale(YScale) ; 

SetXScale (XScale) ; 


if Data = 1 then 

begin 
SetLineStyle(DashedLn, 0,NormWidth) ; 
SetColor (Red); 
DrawLine(YScale, XScale, SetMuF, SetEtaF, SetCF); 
SetColor (Green) ; 
DrawLine(YScale, XScale, SetMuMA, SetEtaMA, SetCMA) ; 

end; 

SetLineStyle(SolidLn, 0,NormWidth) ; 

SetColor (Red); 

DrawCurve (YScale, XScale, NOFVec) ; 

Lane (550, l5y7580,15); 

SetColor (Green) ; 

DrawCurve (YScale, XScale, NOMAVec) ; 

ame (550, 257580,25); 
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SetColor (White); 
OutTextXY (400,470, ’*** Press <Enter> Comer ie eee 
readln; 
CloseGraph; 
end; {Procedure GenDataGraph} 


{+++PART5 : GENERATE COST GRAPH ++4++++++++++++++++++++++++++} 
{ These procedures plot the graph of the mean values and two 
standard deviation bounds for cost advantage of upgrade. } 


Procedure GenGraph; 

var GraphDriver, GraphMode : integer; 
YScale, XScale : real; 
YZero : integer; 


LFHEEHHEEHE THF HEE HEF H+ +++ 4+4+4+4+4+4+4+4+4+4+4+4+4+4} 


Procedure DrawAxis; 
Var 5 = St rang: 


begin {Procedure DrawAxis} 
SetColor (White) ; 
OutTextXY(10,10, ‘Cost advantage’); 


Out TextXY(10,20,° of Upgrade == 7; 

Out Text XY(250,10, ‘Assessing uncertainty’); 
Out Text XY(250,20, ‘from time COLE: 
StrtsT sie 

OUETEXEXKY (350,20, S)%, 

Stree, S): 


OQUETEXEXKY (400, 20,5); 

Line (30,30,307430).; 

Line (30, 430,630,430) ; 
end; {Procedure DrawAxis} 


LH FEEEHEE EEE E HEE FETE HHT HEHEFHE ++ +++ +++4++4+4+} 


Procedure SetYScale (YScale : real; YZero : integer); 
var I, Y ,E: integer; 

CostDiff, M, Scale, Length, Units : real; 

S =; Stung, 


begin {Procedure SetYScale} 

CostDiff := Max-Min; 

M e910; 

jy Sle 

While (CostDif££/M) = b0nde 

begin 
M: 
E 

end; 

Scale += M/Z, 


M0: 
E+1; 
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Hength = 0; 


Units := 0; 

a= 0; 

while Length > -Max do 

begin 
Y := round(Length*yYScale) ; 
Line (30, YZero+Y, 630, YZero+y) ; 
if (I mod 2) = 0 then 
begin 


Str qUnits:1:07S): 
OutTextXY (15, YZero+yY,S); 


end; 
Length := Length-Scale; 
Units := Units+0.5; 
ie s= 1+ 
end; 


Scale := M/2; 


Length := Scale; 
mits := -0.5; 
i= 
while Length < -Min do 
begin 
Y := round(Length*YScale) ; 


Line (30, YZero+Y, 630, YZero+y) ; 
if (1 mod 2) = 0 then 
begin 

Ser (Wmats:1:0,S):; 

Out TextXY (15, YZero+yY,S); 


end; 
Length := Length+Scale; 
Units® -= Units-0.5; 
aa 1: 
end; 
Seats) > 


OutText XY (30,470, 'Y-Axis scale E+’); 
Out Textxy (150,470,S); 
end; {Procedure SetYScale} 


LEEEEEEEEEE EE EEEE HEHEHE EEE EHH t+ t+ t+4+4+4++4+4} 
Procedure SetXScale (XScale : real); 
vom, Tlyex = integer; 


S : string; 


begin {Procedure SetXScale} 


mw -= ST; 

mies = ST mod 5; 

af 1 = 0 then 
mee = T+5—Ri > 


while T <= ET do 
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begin 
X := round ((2=ST) *xSeaq le) ; 
Line (30+X,30,30+X,430); 
if (T mod 10) = O then 
begin 
Str (CES. 
Out Text XY (15+X,440,S); 
end; 
T= Eee 
end; 
Out Text XY (300,460, ‘Month’ ) ; 
end; {Procedure SetXScale} 


LHEEEEEEE EEE EE HEHEE EEE EEE H+ 44¢4444¢4+4+4+444} 
Procedure DrawCurve (YScale, XScale : real); 
var X1, X2, T : integer; 

YM1, YM2, YUB1, YUB2, YLB1, YLB2 : integer; 


begin {Procedure DrawCurve} 


rh: =: Sir. 
Me += 30> 
YM1 := YZero-round (Mean[T]*YScale) ; 
YUB1 ss= YZero-round(UB([T]) *YSeale) ; 
YLB1 := YZero-round(LB[T]*YScale) ; 
for T := (ST+1) to ET do 
begin 
X2 := round( (T-ST) *Xscale+30) ; 
YM2 := YZero-round(Mean[T]*YScale) ; 
YUB2 := YZero-round(UB(T] *YScale) ; 
YLB2 := YZero-round(LB[T] *YScale) ; 


SetColor (Green) ; 

Line (X1,YM1,X2,YM2); 
OutbTextxy (X2=—379M2-5 3a 
SetColor (Red) ; 

lane (X1, YWB1, X2, YUB2) 
Line(X1, YEUBIOX2 | YEe2Z 

Out Text XY (X2-3, YUB2=5, '°.'); 
Out Text XY (X2-3, YLB2-5,’.’); 


ll Ses Ge = 
Vive = YMi2e 
YUB1 := YUB2; 
YEet = YEBZ2- 
end; 
end; {Procedure DrawCurve} 


{HEF HEEEEHEE EEF EHEE HEEFT +++ ++ +44+4+4+4+4+4+4+4} 


begin {Procedure GenGraph} {Verify Driver, Mode, and path} 
GraphDriver := VGA; {adjust to the acceptable values} 
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GraphMode := VGAHi; 
InitGraph (GraphDriver,GraphMode, ‘a:\bgi\’); 
Set BKColor (Black) ; 


ClearDevice; 

DrawAxis; 

YScale := 400/(Max-Min) ; 

XScale := 600/(ET-ST); 

YZero = round(abs (Max) *YScale) +30; 


Line (30, YZero, 630, YZero) ; 
SetLineStyle(DottedLn, 0,NormWidth) ; 
SetyScale(YScale, YZero) ; 
SetXScale(XScale) ; 


SetLineStyle(SolidLn,0,NormWidth) ; 
DrawCurve (YScale, XScale) ; 
faumne (600,25,630,25): 
SetColor (Green) ; 
fine (600,15,630,15); 
SetColor (White) ; 
Out Text XY (500,10, ‘Mean’); 
Out Text XY (500,20, ‘Mean+(-)2SD’); 
Out Text XY (400,470,'*** Press <Enter> to exit ***‘); 
readln; 
CloseGraph; 
end; {Procedure GenGraph} 


{+++PART6 : TWO METHODS & MAIN PROGRAM+4+++4+4+4+4+4+++++4+++++4+++++4+4} 
{ It includes: 

BuildDataArray, 

Likelihood, 

Bayesian, 

Main Program. } 


Procedure BuildDataArray; 


{ This procedure build the array of number of failures 
(maintenance actions) } 


Procedure WriteRealVec (K, H : integer; 
XVec, YVec : RealVec; 
var Output sstext x; 
var T : integer; 


MF : real; 

begin 
writeln(Output, ‘’"Mean number of" "Mean number of" ‘); 
writeln (@utput, ’ Failures Maintenance Actions’); 


writeln(Output) ; 
for T #= K to Hade 
Watt elmOQutput , XVec(T]:12:3,YVec(T] :18:3); 
writeln(Output) ; 
ME «= NS*UPMs 


Jealak 


if Data = 1 then 
for T := H downtowmkK wdeée 
writeln(Output2,MF/XVec[T],MF/YVec[T] ); 
end; 


LFF HEHEHE HEHEHE HEH H ++ +++ 4++4+4+4+4+4+4+4+4+4+4} 


begin {Procedure BuildDataArray} 
Max <:= 90; 
Min := 1E10; 
1f Data = 0 then 
ReadData (Input, NOFVec, NOMAVec) 
else 
begin 
BuildVec1 (ET, SetCF,SetSigmaSaqrF, 
SetMuF, SetEtaF,NOFVec) ; 
BuildVecl1 (ET, SetCMA, SetSigmaSqrMA, 
SetMuMA, SetEtaMA, NOMAVec) ; 
end; 
GenDataGraph; 
WriteRealVec(1,ET,NOFVec, NOMAVec, Output) ; 
end; {Procedure BuildDataArray} 


LFF FEFEEE EE HEEEFEEEEEEEE HEF HE EE EEEE FEE HEE HEH ++ F444 44444444444} 


Procedure Title; 


begin 
if Meth = 0 then 
begin 
writeln(’Maximum likelihood procedure :’); 
writeln(Output, ‘Maximum likelihood procedure :’); 
end 
else 
begin 
writeln(’Bayesian procedure :’); 
writeln (Output, ‘Bayesian procedure :’); 
end; 


writeln(’The best time for subsystem upgrade’ ); 
writeln(’and assessing uncertainty’); 
writeln(’ (frem time ’,ST," to time sr, ja 


writeln; 

write (‘"Time" "Best “ “Cost. “Std 967," 70 
writeln(’ "Mean-2SD" "Mean+2SD"’); 

write (’ Index upgrade advantage Of “COST ).; 
writeln(’ bound bound me 

writeln(’ time of upgrade’); 


writeln(Output, ’The best time for subsystem upgrade’); 
writeln(Output,’and assessing uncertainty’); 


Was 


writeln(Output, 


writeln(Output) ; 


write (Output, 
writeln(Output, 
write (Output, 


writeln(Output, 
writeln(Output, 
end; 


Secommrcime. ~,5T,° to time ’',ET,’)"); 
meme “Best “ " Cost “ "StdeDeve”’); 
: "Mean-2SD" "Mean+2SD"'); 

‘ Index upgrade advantage OL Gest '); 
; bound bound eo); 

t time of upgrade’); 


LA FEAEEEEEEEEEEE FEE FEE EEE EEEFE FEE EEF H FEF EHH HHH 444444444444} 


{ This procedure applies maximum likelihood model. } 
Procedure Likelihood; 
begin 

Max := 0; 

Min := 0; 

Title; 

for Time := ST to ET do 


begin 


FindHat Parameter1 (Time, HatCF,NOFVec,MuF,EtaF, 


SigmaSqrF,HatMuF, HatEtaF) ; 


FindHat Parameter1 (Time, HatCMA, NOMAVec,MuMA, EtaMA, 


SigmaSqrMA, HatMuMA, HatEtaMA) ; 


FindChangeTimel (Hor, Time+LeadTime, HatCF, HatCMA, 


Total := 0; 


HatMuF, HatEtaF, HatMuMA, HatEtaMA) ; 


Reset Parameterl1; 
if Rep > 0 then 


— 


tor R 
begin 


1 to Rep do 


BuildVec (Time, SetCF, SetSigmaSarF, 


SetMuF, SetEtaF, SimuFVec) ; 


BuildVec (Time, SetCMA, SetSigmaSqrMA, 


SetMuMA, SetEtaMA, SimuMAVec) ; 


FindHat Parameter1 (Time, HatCF, SimuFVec,MuF,EtaF, 


SigmaSqrF, HatMuF,HatEtaF) ; 


FindHat Parameter1 (Time, HatCMA, SimuMAVec,MuMA, 


EtaMA, SigmaSqrMA, HatMuMA, 
HatEtaMA) ; 


Uncertaintyl1 (Hor, Time+LeadTime,R,HatCF,HatCMA, 


end 
else 
begin 
writeln; 


HatMuF, HatEtaF,HatMuMA, HatEtaMA) ; 


writeln (Output) ; 


end; 
end; 
writeln (Output) 


° 
f 


LS) 


end; 


LAFAFEFEE FEF EEE EE FEF EEE EEE EEE EE FE FE FEF FEE EHH FH 4 44444444444} 
{ This procedure applies the Bayesian model. } 


Procedure Bayesian; 
{H+ t+ ++ +t +++ ++ t+tt+ttett+t+tt+++ttt+++t+} 


Procedure Draw(var HatC : integer; 
var HatMu : real; 
var HatEta : real; 
Mu, Eta, RSqr, VSqr, RhoSqr, PiDist : RealVec); 
var C : integer; 
U, =f, COV 9 real 


begin 
U s= Random; 
Cae ls 
whilest) > PiDist[C] do 
C :=3er1- 
HatCee= Ge 
Z := GenNormal; 


HatMu := Mu[C]+(sqrt(RSqr[C])*Z) ; 
if C =’ Time then 


begin 
Cov := (-=S@met(Raesaqn |e) 74 
+(sqrt (1-RhoSqr[C] ) *GenNormal) ; 
Hatta := Bta[C]+(sqrt (VSqr([C] ) *@ovwer 
end 
else 
HatEta :="0; 
end; 


LEE EHEEE HEHEHE HEHEHE +H ++ ++ t++++e++++++i 


begin 
Max 
Min 
Ti the: 
InitialSigmaSqr (NOFVec,MuF,EtaF,SigmaSqrF) ; 
InitialSigmaSqr (NOMAVec, MuMA, EtaMA, SigmaSqrMaA) ; 
fer Tame :=SSm to ET doe 
begin 
FindHat Parameter2 (Time, NOFVec,MuF, EtaF,SigmaSarF, 
VSqrF,RSqrF, RhoSqrF, PiStarF) ; 
FindHat Parameter2 (Time, NOMAVec,MuMA, EtaMA, SigmaSqrMA, 
VSqrMA, RSqrMA, RnoSqrMA, PiStarMaA) ; 
FindChangeTime2 (Hor, Time+LeadTime, MuF, EtaF, 
MuMA, EtaMA, PiStarF, PiStarMA) ; 


OW 
Oi; 
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Lota = 0; 
Reset Parameter2; 
if Rep > 0 then 
Uncertainty2 (Hor, Time+LeadTime, MuF, EtaF,MuMA, EtaMA, 
RSqrF, VSqrF,RhoSqrF,RSqrMA, VSqrMA, 
RhoSqrMA, PiStarF, PiStarMA) 
{ For R 3= 1 to Rep do 
begin 
Draw (HatCF,HatMuF, HatEtaF,MuF, EtaF, 
RSqrF,VSqrF,RhoSaqrF, PiDistF) ; 
Draw (HatCMA, HatMuMA, HatEtaMA, MuMA,EtaMA, 
RSqrMA, VSqrMA, RhoSqrMA, PiDistMA) ; 
Uncertaintyl (Hor, Time+LeadTime,R,HatCF,HatCMA, 
HatMuF, HatEtaF, HatMuMA, HatEtaMA) ; 


end} {You may use the procedure in the pair of } 

else {braces instead of using procedure } 

begin {Uncertainty2. The procedure in the braces} 

writeln; {uses the resampling method to assess } 

writeln(Output); f{uncertainty. } 

end; {Be sure to put a pair of braces on one of} 

end; {the two procedures. } 
writeln(Output) ; 


end; 
LHEEEEEEEEFR HEE EEE EEE EEEEEHEEEEPHEEE EE HEE EEE HEE EE PE +E HH + tet te+44} 


Procedure ViewGraph; 
marmeCh : char; 


begin 
writeln; 
writeln(’Do you want to see the graph’); 
weate(’<y/n>? ‘); 
readin(Ch) ; 
Meme = 1°) Or (Ch = ‘y’) then 
GenGraph; 
end; 


LHEEEEEEEEEEFE EEE EEF EEEEE EE EEEE EE FE FEET EEF HET E HEH +++ 4+4+4+4+4+4+4 } 
{ This procedure gives the choise to reuse the simulated 
data. } 


Procedure Again; 
wameen > char; 


begin 
writeln; 
writeln(’Do you want to use the simulating data again’); 
writeln(’Data will be contained in "Simu.data" file’); 
write(‘’<y/n>? ‘); 
readin(Ch); 


dS 


1f (Ch = °Y') Sem Gh =F weewen 


begin 
reset (Output2) ; 
Data =e): 
Infile := ‘a:\Simu.Data’ 
end; 
end; 


LAFEEEAEEEEEALHEAEEE EEE EEE EEEEHHEEH HEE HEF HE ++ +++ tt+t+t++t+++++++4 i 
{ This Function let user decide to terminate or continue the 
Job. } 


Function Done : boolean; 
var Ch : char; 


begin 
repeat 
writeln; 
writeln(’Enter 0 to stop or 1 to continue.’); 
WIriteuaren®) ; 
readln(Ch) ; 
Mien (eh = 0 ") 
Done -:=s(Ghr= '0’ 
end; 


Oe meee ee 
) O 


c 


LEEEEEEEFEEEEEE EEE EE EEFEFEFEFEFEFEPTHEFE HEHEHE ++ Ht ++ t+ tt++t+++++ i} 


begin {Main Program} 
assign (Output,’a:\F.out’); 
assign (Output2,’a:\Simu.Data’) ; 
rewrite (Output) ; 
rewrite (Output2) ; 


DefaultParameter; 
SetUnitCost; 
SetMission; 


Randomize; 
RanSeed := 1; {User may change the value of RanSeed} 


repeat 
SetParameter; 
WriteParameter; 
BuildDataArray ; 
if Meth = O then 
Likelihood 
else 
Bayesian; 
if Rep <> 0 then 
ViewGraph; 
if Data = 1 then 


aa 


Again; 
until Done; 


close(Output) ; 
close(Output2) ; 
end. {Main Program} 
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